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Abstract. For the wave equation associated to the Laplacian on a compact 
manifold with boundary with a conic metric (with respect to which the bound- 
ary is metricaUy a point) the propagation of singularities through the boundary 
is analyzed. Under appropriate regularity assumptions the diffracted, non- 
direct, wave produced by the boundary is shown to have Sobolev regularity 
greater than the incoming wave. 



Introduction 

Solutions to the wave equation (for the Fricdrichs extension of the Laplacian) 
associated to a conic metric on a compact manifold with boundary exhibit a diffrac- 
tive, or "ringing," effect when singularities strike the boundary. The main results 
of this paper describe the relationship between the strength of the singularities in- 
cident on the boundary and the strength of the diffracted singularities. We first 
show that if no singularities arrive at the boundary at a time t then the solution is 
smooth near the boundary at that time, in the sense that it is locally in the inter- 
section of the domains of all powers of the Laplacian. We then show that if there 
are singularities incident on the boundary at time i and in addition the solution 
satisfies an appropriate nonfocusing condition with respect to the boundary, then 
the strongest singularities leaving the boundary at that time are on the geometric 
continuations of those incoming bicharacteristics which carry singularities, whereas 
on the diffracted, i.e. not geometrically continued, rays the singularities are weaker. 
If the incident wave satisfies a conormality condition, then the singularity on the 
diffracted front is shown to be conormal. Applying this analysis to the forward 
fundamental solution gives an extension of results of Cheeger and Taylor ([2, 3]) 
from the product-conic to the general conic case. The results contained in this 
paper represent a refinement of those previously announced in [15]. 

The problem of diffraction is an old one, with rigorous treatment stretching 
back to the work of Sommcrfeld, who, in 1896, discussed diffraction around edges 
in the plane [22] ; this includes the case of Dirichlet boundary conditions at the slit 
[0, oo) C C. Using the method of images, this problem may be reduced to the study 
of the wave equation on the cone over the circle of circumference 47r. An overview 
of this and many other problems of diffraction around obstacles in R" is given in 
Friedlander [5]. 
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Diffractive effects were extensively studied by Cheeger and Taylor [2, 3] in the 
special case of product- conic metrics, i.e. M+-honiogeneous metrics on (0, oo) x Y 
of the form 

g = dx^ + x'^h{y, dy). 
Separation of variables can be used to give an explicit description of the funda- 
mental solution in terms of functions of the tangential Laplacian. Sec also the 
discussion by Kalka and MenikofF [11]. Roulcux [19] obtained a version of the re- 
sults of Cheeger- Taylor in the analytic category. Lebeau [12, 13] has also obtained 
a diffractive theorem in the setting of manifolds with corners in the analytic cat- 
egory, and Gerard-Lebeau [6] have explicitly analyzed the problem of an analytic 
conormal wave incident on an analytic corner in M? . Estimates in spaces have 
been obtained for product-conic metrics by Miiller-Seeger [18]. 

Let X be an n-dimensional conic manifold, that is to say, a compact manifold 
with boundary, with a Riemannian metric g on the interior which, near the bound- 
ary, takes the degenerate form 

(1.1) g = dx^+x^h, 

where a: is a boundary defining function and h e C°°(X; Sym^(T*X)) restricts to 
be a metric on dX. Each boundary component of X is thus a "cone point" in the 
metric sense. A trivial example of a conic metric is obtained by blowing up a point 
p'm& Riemannian manifold; near p we can take x to be the distance function to p, 
and (I.l) is simply the expression of the metric in Riemannian polar coordinates. 
Such conic metrics exist on any manifold with boundary. 

The subject of this paper is the wave equation on a conic manifold, and in par- 
ticular the propagation of singularities for its solutions. Let A be the Priedrichs 
extension of the (non-negative) Laplace-Beltrami operator on X. We consider so- 
lutions to 

(1.2) {Dl - A)u = 
on R X X° which are admissible in the sense that 

u e C"°°(M; Dom(Zi"/2)) ^j. ^^^^ ^ 

and the equation holds in C-°°(K; T)ora{A''/'^'^)). 

Hormander's theorem [9] on the propagation of singularities for operators of real 
principal type yields rather complete information about the location of the singu- 
larities of any solution of (L2) away from dX. Namely, WF''(m) (the wavefront set 
computed with respect to the scale of Sobolev spaces) is contained in the charac- 
teristic variety and is a union of maximally extended null bicharacteristics. Since 
the null bicharacteristics are essentially time-parametrized geodesies over X° , this 
means, somewhat loosely speaking, that the singularities travel with unit speed 
along geodesies. 

Hormander's theorem does not address the question of what happens when sin- 
gularities reach dX. That is, it leaves open the question of how singularity-carrying 
geodesies in dX terminating at dX are connected with those emanating from dX . 
The answer involves both the geometric propagation of singularities that one would 
expect from the limiting behavior of geodesies that come close to the tip of the cone 
without striking it, and an additional, diffractive, effect. 

By the finiteness of the propagation speed for the wave equation it suffices to 
consider a single component Y of dX (i.e. a single metric cone point). We first 
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note that for every y G Y, there is a unique (maximal, unparametrized) geodesic 
in X° — X \ dX with y in its closure; conversely there is a neighborhood of Y in 
which each point is the end point of a unique, short, geodesic segment terminating 
at Y. This corresponds to a product decomposition [0, e)^ x y of a neighborhood of 
Y in which projection to Y is the map to the end-point and x is the length of this 
segment. We henceforth take x to be the boundary defining function; the metric 
still has the form (I.l) but now with h{dx, ■) = 0, so /i is a family of metrics on Y 
parameterized by x. 

It follows that the surfaces t ±x = c, for any constant c, which are well-defined 
near each boundary component, are characteristic for the wave equation. These 
radial surfaces carry the extra, diffractive, singularities, the existence and regularity 
of which is our main object of study. 

Our results are best illustrated by the fundamental solution itself. Consider an 
initial point m G X°, which we shall take to be close to some boundary component 
Y. There is then a unique short geodesic interval from fh to the boundary; its 
length is rf(m,F). The wave cone emanating from m is smooth for small positive 
times but for longer times generally becomes singular. However it is the projection 
of a smooth Lagrangian submanifold Im C T*(R x X°) consisting of the union of 
those maximally extended null bicharacteristics of the wave operator which pass 
above m at t = 0. In addition we consider the radial surface mentioned above, 
Dfh = {t = X + d{fh,Y)}, which is well defined if d{fri,Y) is small enough and 
which emanates from the boundary at the time of arrival of the geodesic from fh. 

Theorem I.l (Fundamental solution). Let Efn be the fundamental solution of the 
wave equation for the Friedrichs extension of the Laplacian of a conic metric on a 
compact manifold with boundary X, with pole at fh € X°. If d(fh, Y) is sufficiently 
small then, regarding as a distribution on M x X°, 

(1.3) WF(i;^) C 7^ U N*Drn in \t\ < 2d(m, Y) 

and Efn is conormal with respect to Dm away from the projection of the closure of 
Ifji and, is of Sobolev order \ — ^ there for any (5 > 0. 

A comparison with the explicit results of Cheeger- Taylor [2, 3] in the special case 
of product cones shows that this result is optimal as far as Sobolev regularity is 
concerned. 

The first step in proving this is to show that no singularities arise from the 
boundary spontaneously. 

Theorem 1.2 (Diffractive Theorem). // an admissible solution to the wave equa- 
tion on a conic manifold has singularities, of Sobolev order r, on at least one null 
bicharacteristic hitting a boundary component at time i then it has singularities, of 
Sobolev order r, on at least one null bicharacteristic leaving that boundary compo- 
nent at time t. 

This diffractive theorem does not give any localization of singularities in the bound- 
ary component Y. It in no way distinguishes between the different "outgoing" null 
bicharacteristics leaving the boundary component at t = i although in fact it is 
easily strengthened to show that singularities in the two components of the char- 
acteristic set, corresponding to the sign of the dual variable to t, do not interact at 
all. 

In order to state a more refined theorem, which distinguishes between different 
points in Y, we need first to consider the possible geometric continuations of a 
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geodesic terminating at y e F. Geometric continuation of geodesies corresponds to 
the relation 

(1.4) 

G{y) — {y' Cz Y;y and y' are endpoints of a geodesic segment of length tt in Y}, 

where Y is endowed with the metric h\x=o ■ The geodesic segments of length tt in 

Y arise naturally as limits of geodesies in X° which narrowly miss Y. 

Our second main result shows that under certain circumstances, a singularity 
arriving at a point y GY a,t t = t will produce only weaker singularities along rays 
emanating from points in y \ G{y) than on the geometrically continued rays, which 
are those emanating from G{y). In other words, we show that diffracted singularities 
are weaker than geometrically propagated ones. 

The further hypothesis that we make is one of nonfocusing of singularities. Sup- 
pose that the solution u is microlocally of Sobolev order r on all incoming null 
bicharacteristics arriving at Y at t = i. We require that tangential smoothing im- 
proves this regularity; for our purposes the amount of tangential smoothing required 
to gain regularity is irrelevant. Let Ay denote the Laplace-Beltrami operator on 

Y with respect to the metric h\x=o, extended to operate on a collar neighborhood 
of Y by the metric product decomposition described above. We assume that for N 
sufficiently large (1 -|- Ay)~^u is more regular (microlocally) than u in the sense 
that 

(1.5) p^WF''' {{1 + Ay)~^u) iorpon 

the incoming null bicharacteristic reaching Y a,t y gY, t = t. 

Theorem 1.3 (Geometric Propagation). Ifu is an admissible solution of the wave 
equation, then the outgoing null bicharacteristic emanating from a point y & Y 
at t = i is disjoint from WF^(u) provided all the incoming null bicharacteristics 
reaching G{y) at t = t are outside WF^('u) and (1.5) holds with r' > R microlocally 
for all other incoming null bicharacteristics meeting the boundary at time i. 

This theorem may not be strengthened by dropping the nonfocusing assumption 
(1.5) as is shown by a counterexample in §16. On the other hand the two components 
of the characteristic variety are again completely independent. 

At a time i the background regularity, r, for an admissible solution near a bound- 
ary component is the Sobolev regularity that holds on all incoming null bicharac- 
teristics arriving at that boundary component at t = i. Theorem 1.2 shows that 
this background regularity propagates to all outgoing null bicharacteristics leaving 
the boundary at time t. Theorem 1.3 above shows that there is additionally a gain 
of any number, I, of Sobolev derivatives on any one of these outgoing bicharac- 
teristics, say emanating from y G Y, provided that u has Sobolev regularity r + I 
on all incoming bicharacteristics which are geometrically related to this one and 
that u has the property that sufficient tangential smoothing increases its regularity 
microlocally near all other (and hence all) incoming null bicharacteristics, reaching 
the boundary at time i, to greater than r + 1. 

Although there is in principle no upper limit to the gain of regularity compared 
to backgroimd available from Theorem 1.3, there are practical limits. The most 
important case of nontrivial tangential smoothing is when the solution is, for t <t 
and away from the boundary, a conormal distribution associated to a hypersurface 
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that is simply tangent to the incoming radial surface x = t — t. For such a distribu- 
tion, the lemma of stationary phase shows that tangential smoothing gains almost 
(n — l)/2 derivatives, hence Theorem 1.3 guarantees that the diffracted wave is 
almost (n — l)/2 derivatives smoother than the main singularity. In addition, the 
diffracted wave for such a solution will itself be conormal. 

More precise statements of the theorems above as well as further results are given 
in Section 4. We now briefly describe the ingredients in the proofs. 

Wc make extensive use of the calculus of edge pseudodifFerential operators as de- 
veloped by Mazzeo [14]. A related calculus of pseudodifferential operators adapted 
to edge structures has also been constructed by Schulze [20] but the lack of "com- 
pleteness" in this calculus makes it much less applicable; in particular it does not 
seem to have associated with it a useful notion of wavefront set. We develop and 
use just such a notion for Mazzeo's calculus measuring microlocal regularity with 
respect to the intrinsic weighted Sobolev spaces. We use this to obtain a result on 
the propagation of singularities at dX analogous to Hormander's interior result. 
The technique of the proof is a positive commutator construction similar to that 
used in [9], but with the crucial distinction, familiar from scattering theory, that 
the bicharacteristic flow now has radial points, thus necessitating a more subtle 
construction. At the radial surfaces the propagation arguments are only valid for a 
limited range of weighted (edge) Sobolev spaces, corresponding to a linear estimate 
on the weight in terms of the regularity. Such an estimate amounts to a divisibility 
property for the solution in terms of powers of the boundary defining function. 

The proof of the diffractive theorem. Theorem 1.2, relies on the extraction of 
a leading part, the non-commutative normal operator in the edge calculus, via an 
appropriately rescaled FBI (Fourier-Bros-Iagolnitzer) transform S near Y. The 
model operator is Z\o ^ l, where Aq is the Laplacian on the tangent cone R+ x Y 
for the product-conic metric dx^ + x^{h\x=o)- Iterative application of the outgoing 
resolvent for this model operator, corresponding to the scattering structure at the 
infinite volume end of the cone, combined with the microlocal propagation results 
discussed above, yields the regularity of u. We use estimates from [16], although in 
this product case we could instead rely on direct methods, as used for instance in 
[2]. 

To prove the geometric propagation theorem, we begin by establishing a division 
theorem to the effect that the nonfocusing assumption (1.5) implies that the solution 
has better decay in x than would be predicted merely by energy conservation. This 
allows us to apply the propagation results to prove a special case of Theorem 1.3. 

We also show that conormality on a radial surfaces persists for solutions. That is, 
if (w(0), Dtu{0)) are conormal distributions at a surface x = x, supported away from 
a boundary component Y , and if x is sufficiently small then by standard interior 
regularity results {u{t) , Dtu{t)) , is conormal with respect to {x = x — t} D {x = 
X + t} for small positive t < x; we show that it continues to be conormal with 
respect to {a; = \x — 1\} U {x = x + t} for small t > x, i.e. after the wave leaves 
the boundary. Such conormal solutions provide the counterexamples mentioned 
following Theorem 1.3. 

For solutions, such as the fundamental solution, that are initially conormal at 
a surface meeting the radial surfaces x = x — t with at most simple tangency, 
we are additionally able to prove that the diffracted front is conormal; this is in 
effect a microlocal version of the radial conormality argument discussed above. The 
conormality of the diffracted front and the special case of Theorem 1.3 mentioned 
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above suffice to establish conormality and (sharp) regularity for the diffracted wave 
of the fundamental solution, which is then used to prove Theorem 1.3 in its full 
generality. 

The outline of the paper is as follows. In §1 we prove the existence of a product 
decomposition for a conic metric in a neighborhood of dX; this is equivalent to 
reducing the metric to the normal form 

g — dx'^ + x^h{x, y, dy). 

Then in §2 and §3, we discuss the mapping properties and domains of the Laplacian 
and its powers on a conic manifold. This enables us to give precise statements of the 
main theorems in §4. Sections 5-8 build up the machinery of the edge calculus and 
culminate in the proof of the propagation theorem for the edge wavefront set. Then 
in §9 and §10 we discuss the rescaled FBI transform and the normal operator in the 
edge calculus, which wc use to prove Theorem 1.2 in §11. In §12 we demonstrate 
the conservation of tangential regularity, i.e. iterated regularity under powers of 
Ay- This is then used in two ways. In §13 we use conservation of tangential 
regularity together with conservation of iterated regularity under vector fields of 
the form (xD^+tDf) to prove conservation of radial conormality. In §14, we use the 
conservation of tangential regularity to prove the division theorem; the proof of the 
sharper division result for conormal initial data is closely related to that of radial 
conormality. In §15 we prove a preliminary version of Theorem 1.3 and establish a 
theorem on conormality of the diffracted front. In §16 we discuss a class of examples 
including the fundamental solution and prove Theorem I.l; as a consequence we 
then prove Theorem 1.3 in full generality, and discuss a counterexample to that 
theorem when the nonfocusing condition is omitted. 

The authors thank Daniel Grieser and Andras Vasy for comments on the man- 
uscript. 

1. Conic metrics; geodesics and normal form 

Let X be an n-dimcnsional manifold with compact boundary, dX its boundary 
and X° = X \ dX its interior. 

Definition 1.1. A conic metric on X is a Riemannian metric g on X° such that 
in a neighborhood of any boundary component Y oi X, there exists a boundary 
defining function x {x > 0, {x = 0} = dX, dx\Qx^ 0) in terms of which 

(1.1) g = dx^+x^h 

where h G C°°(Sym^T*X) and ft, fy is a metric. A conic manifold is a compact 
manifold with boundary endowed with a conic metric. 

Let F be a compact manifold without boundary. The product- conic metrics on 
[0, oo)x X Y are the conic metrics of the form 

(1.2) g = dx^ + x^ho 

where /iq is a metric on Y. These are the model cases which form the basis of our 
analysis below. They also motivate the basic normal form for conic metrics which 

we discuss next. 

Theorem 1.2. Let g he a conic metric on X. There exists a collar neighborhood 
O of dX and an isomorphism {x, 11) : O — > [0, e) x dX, e > 0, such that in terms 
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of this product decomposition 

(1.3) g = dx^+ x^JThix) near dX, 
where h{x) G C°°{[0, e); metrics on dX). 

Proof. The existence of such a normal form distinguishes a vector field V by 

(1.4) Vx = l and Lvh{x) = near dX. 

The integral curves of this vector field are geodesies. Our main task is thus to show 
that through any point in a neighborhood of dX, there is a unique short geodesic 
reaching dX\ the length of this short geodesic segment will then furnish the desired 
defining function for dX. 

As our construction is local near a boundary component, we assume without loss 
of generality that dX is connected. We begin with any collar neighborhood and 
associated projection -kqx ■ O dX, and choose a boundary defining function p 
for dX and coordinates v = ttqxV with y any coordinates on dX. We may asstune 
p to have been chosen in accordance with our definition of a conic metric, so that 

g = dp^ + p^h 

where /i is a smooth, symmetric two-tensor. Let 

h = •dp^ + •dpdv + kij{v)dv^dv^ + P* dv^dv^; 

then 

(1.5) (•, = (1 + o{p^))dl + o{i)dpd^ + {p-^¥i{v) + o{p-^))d^,d^, , 

where, here and henceforth, 0{p^) means p^ times a smooth function of {p,v). 

We now regard the dual metric as a function on a rescaled version of the cotan- 
gent bundle. Let ^TX denote the b-tangent bundle whose sections are the vector 
fields tangent to dX. Let '°T*X be its dual. Then sections of ^T*X are C°°-linear 
combinations of dp/p and duj's. Writing the canonical one-form on ^T*X as 

(1.6) ^^+r)-dv, 

P 

we can now write the dual metric as 

(1-7) ((^, ry), (^, nX = ^-t^^ + 0{l)e + Oip-')^v + 0{p-'W- 

" p 

Over the interior '^T*X° = T*X° and the canonical symplectic form on T*X 
lifts to a form on ^T*X, singular at the boundary, given in canonical coordinates 
by d{£,dp/p + rj ■ dv). Associated with this symplectic form and the energy function 
(•, •)* is the Hamilton vector field 

Hg = ^i^pdp + ie + k'^niVj)di + Hex) + P, 



where 



1 

7 



P = -{{0{p^)i + 0{p')7^)dp + {0{p^)e + 0{p)£,7i + 0{p)Tf)d^ 



+ {0{p)i + 0{p)^)d^ + {0{p')e + 0(p)er; + 0(p)r;^)5„ 

and Hqx is geodesic spray in the {v, rj) variables with respect to the metric k. The 
projections of the integral curves of Hg to X are geodesies. 
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Let Hq = Hg — P he the main term in the Hamilton vector field. The vector 

field _ffo is thus tangent to the submanifold TV = {77 = 0} of "normal" directions to 
dX. Since on N, the flowout of {p^ /2)Hq \s p = po/(l — ios) and ^ = ^o/(l — £,os), 
the map 

p I— > lim exp 

s— »+oo 

maps to dX smoothly in a neighborhood of the boundary, takings 1— > v{p) e dX. 
The projections of integral curves of i?o on N would be geodesies reaching the cone 
point were the metric of the desired form (1.3). 

We now show that an analogous normal manifold, given by a perturbation of N , 
exists for the full vector field Hg. The flow along this manifold to the boundary 
will provide the desired geodesies to the cone point. 

To simplify the discussion of the vector field, Hg, we scale away its homogeneity 
in vi) as follows. Consider the smooth function 

1 

on the complement of the zero-section in ^T*X. Note that <j is approximately equal 
to (^^ + k^^Tiirfj)'^/"^ , and is homogeneous of degree —1 in the fibers of ^T*X. We 

further set 

(1.8) f = 77 = ^V- 

Since (•, •)* is preserved under the flow Hg, {x^ /2)Hga = + O^p"^)^ + 0{p)fi. 
Note that the error terms in this expression are at least quadratic in {p,fj). Thus 
we can write 

(1.9) 

1 - - ■• •• /- I dk^^ \ 

-p'^aHg = ^ada + ^pdp + k'^f]jd^, + k'^fjifjjd^ - \ ik + Ydv^^'^A ^ ^' 

where the perturbation term has the form 

P' = 0{p)id,, + 0{p^ + f) (C°° vector field) ; 

the single non-quadratic error term comes directly from the corresponding term in 
P. By homogeneity, the error terms above are independent of a, i.e. {l/2)p'^aHg 

pushes forward to a vector field on the unit sphere bimdle ^S*X = ^T*X/M.+. On 
^S*X, {p,v,fi) are coordinates near f] = 0, since 

{{U),{^,v))*g = 1 at 77 = 0, 
^ = ±1. Moreover we easily see that ^ T 1 vanishes to second order at p = fj = 0; 
of crucial importance is the fact that the power of p multiplying the error term 
in (1.7) is larger than those multiplying the ^rj and 77^ terms. 

In these coordinates, the linearization of the vector field {l/2)p-^aHg in {p,f]) 
near p = fj = and ^ = ±1 is simply 

L = ±pdp + {k'^fjj + 0{p))d^, T 

Since {l/2)p'^aHg vanishes identically at p = 77 = 0, and its linearization has 
eigenvalues ±1 in the normal directions, {1/2) p^aHg is r-normally hyperbolic near 
p = 7^ = for all r, in the notation of [7]. Hence by the Stable/Unstable Manifold 
Theorem as stated in Theorem 4.1 of [7], near {p = = 0, |^ = 1} there exists a 
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stable invariant manifold N for the flow of {l/2)p'^aHg, with N tangent to {t] = 0} 
at {p = fj = 0}. Because N is tangent to = 0}, the projection map 

TT : ^S*X X 

restricts to give a diffeomorphism N = X (locally, for p sufficiently small) . 

On N, lims_>_ooP = 0, where s parametrizes the flow along {l/2)p'^aHg; using 
p as a parameter along the flow on A'' yields, by (1.8) and (1.9), 

- = o{i) + oifj/py, 

the latter term is in fact 0(1) since N is tangent to {fj = 0} at p = 0. Thus, 
lims^_oo V is in fact a smooth map. Hence 

(1.10) n : p ^ hni^exp^-i(p) ^p'^aHg 

is a smooth map from a neighborhood O' of dX in X to dX. Indeed 11 is a fibration, 
with fibers R+ given by the projections of integral curves of {l/2)p^aHg, i.e. by 
geodesies hitting the "cone point" dX. 

On O' , set yi{p) = Vi{Il{p)) and x{p) = dg{p,'n.{p)). To first order at dX, x = p 
and y = v. Hence {x,yi, . . . ,?/„_i) form a coordinate system on a neighborhood 
O". In these coordinates, {dx,dx)g = 1 since exptdx is unit speed flow along 
geodesies reaching dX. Furthermore, {dx,dy.)g = for all i by Gauss's Lemma^ 
and (dy^ , dy^ )^ = 0{x'^). Hence in the coordinates (a;, y), g takes the form (1.3). □ 

Since x is the distance along the normal geodesies it is uniquely determined 

by (1-3); the vector field V determining the product decomposition is also fixed 
geometrically by (1.4). The choice of a conic metric additionally induces a metric 
on dX, namely 

ho = h{0). 

Henceforth, x will always denote this distance function for the given conic metric. 

The proof of Theorem 1.2 used in a crucial way the existence of a unique normal 
geodesic starting at each point of the boundary. These geodesies foliate a neighbor- 
hood of the boundary, and indeed, the existence of this foliation characterizes conic 
metrics among more general nondegenerate forms in dx and xdyi. On IR x X°, 
we thus obtain a foliation of a neighborhood of M x dX by projections of null 
bicharacteristics for the symbol of the d'Alembertian, — A. 

We shall fix notation for various sets corresponding to the geodesic segments 
that hit the boundary and are within small distance e > of it. 

Definition 1.3. Let = {adt + pdx : \a\ = |/3|, a; < e} C T*{R x X°). For ieR, 
yeY, let 

R% j{i,y) = N^n {sgna/J = 1, sgna = ±, t = i-x, y = y}, 
R% o{i^y) = n {sgna/J = 1, sgna = ±, t = t + x, y = y). 

^Thc usual proof of Gauss's Lemma using the first variation fommla works even at a cone point, 
i.e. geodesies with one endpoint on dX are orthogonal to the hypcrsuifaccs dg(p,dX) = c. We 
simply let yiu) be a curve in dX, and apply the first variation formula to the family of geodesies 
connecting x = xq, y = y{u) to x = 0, y = y{u). 
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If any of the parameters • = 7 or O, ±, t, or y is omitted, the resulting set is defined 
as the union over all possible values of that parameter. If F is a single component 
of dX, let 

veY 

The set N'^ is the normal set near the boundary of X of points in the cotangent 
bundle which lie along geodesies (projections of bicharacteristics) entering and leav- 
ing dX. Here "I" and "O" stand for the "incoming" and "outgoing" components 
of N, on which dx/dt is respectively negative and positive. The additional sign 
± microlocalizes in the sign of the dual variable to t. The sets i?!^ ,{t,y) are the 
points through which the short geodesic to the boundary arrives at or departs from 
dX at time t = t and at the point y. 

If X is a compact conic manifold then every geodesic starting at an interior 
point can be extended maximally in both directions until and unless it terminates 
at the boundary. This naturally suggests the question of how such geodesies can, 
or should, be further extended. 

Definition 1.4. By a limiting geodesic in a conic manifold we mean a continuous 
piecewise smooth curve c : I — > X, where / = Ujlj is decomposed as a locally 
finite union of relatively closed subintervals on each of which c restricts to be a 
smooth curve cj and such that 

1. Each Cj is either a geodesic in X or a geodesic (for ho) in dX and such 
segments alternate. 

2. Boundary segments are of length at most tt and if such a boundary segment 
is not the first or last segment then its length is exactly tt. 

Lemma 1.5. If is a sequence of geodesies in X° which converges uniformly as 
curves in X then its limit is a limiting geodesic and conversely an open neighborhood 
of each boundary segment of any limiting geodesic arises as such a limit. 

This result will not be used except as motivation, hence we omit its proof. 

In view of this behavior of the geodesies we define a singular relation, interpreted 
as a set-valued map 

cRqX R} 

such that 

(1.11) for p e i?i,o©, r^(p) = {9 e R^^it) u{q) e G{u{p))}, 

where 11 is the projection map from a collar neighborhood of dX to dX given by 
(1.10), and for y G dX, G{y) is the geodesic relation defined in (1.4). Clearly, 

(1.12) i?l,/(i)= U 

pei?.±,o(t) 

In fact T'^(p) is generically of codimension one in j{t). 

To the given conic metric and a boundary component Y we associated the lim- 
iting product metric 

(1.13) go = dx^+x^ho 

on the normal bundle to Y, which we may identify with [0,oo) x Y using the 
decomposition associated with (1.3). Let A denote the (nonnegative) Laplace- 
Beltrami operator with respect to the metric g and Aq that with respect to go. 
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Then, letting {wq, Wn-i) = {x,yi,..., Vn-i), 

n 

j,k=0 

becomes, near Y, 



2 i[{n-l) + xe] i Ah 



(1.14) A = Di-^ ^_^D,--D, + -^, 



X 



laiogdet/ifa;) 1 /, , , ,dh(x) 
where e = ^ — = 7:tr( h~^ix) ^ ' 



2 dx 2 \ ^ ' dx 

Here Ah is the Laplacian on dX with respect to the (a;-dependent) metric h{x) on 
dX. Similarly 

(1.15) A, = Dl-':^D. + ^, 

x x^ 

hence if we identify [0, oo) x Y with the metric product decomposition near y, 

(1.16) - Zio e x-^ (C°°(a;, y){xD,) + ^ C°°(a;, v)Dy,Dy^ + ^ C°°(a;, y)£>j,,) 

We will also use the notation Ay for Ah^ , when restricting our attention to a single 
boundary component. 

If 5 is a product-conic metric, it is easy to check that 

[□.^y] -0, 
[□, [xD, + (t - t)Dt)\ - ~2iU. 

In the general conic case, these "symmetries" are broken. It is crucial for our 
purposes, though, that perturbed versions of the above identities still hold. In 
particular, if for brevity we set 

(1.17) R = xD^^{t-f)Du 
then 

(1.18) [□, Ay\ = QD^ + x-^P, with 

Q e C°°([0,e);Diffi(y)), and P e C°°([0, e); Diff3(y)) 

and 

(1.19) p,^] = -2in + aL»^ + a;"ip,with 

aGC°°([0,e) X Y) and P G C°°([0, e); Diff2(y)). 

More generally, 

Lemma 1.6. For any q,k gN, 

(1.20) UB!" = J2 Ck,jR^O + J2 (akjD^ + -Pkj J R^, where 

Ckj e C, ak,j e C°°([0, e) x Y) and Pkj e C°°([0, e); Diff2(y)), 
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(1.21) DA'^ = A-^n + Y.iQi^rD, + ^^9,r)^y 

With Qq,r e C°°([0, e); Diff«-''(y)), Pg,^ e C°°([0, e); Diff«-''+2(y)), 



anrf 



(1.22) □Zi«,i?'= =^CfcjZi«.i?^n+ iQk,l,<,,rD:, + ^Pk,l,g,r)^VR' 

i=0 0<r<(j, 0<j<fc 

r+j<k+q 

where Quj,,,r e C-([0, e); Diff^-'-(r)), Pfe,,,,, e C-([0, e); DifF^-'-+2(r)). 

Furthermore, all the differential operators Pkj, Pq,r, Pk,i,q,r have vanishing con- 
stant terms. 

Proof. For k = 1 (1.20) follows directly from (1.19). For general k we may use this 
together with (1.20), as an inductive hypothesis, to see that 

□P'^+i = RUR'' - [R, □jP'^ 

is of the stated form. Equation (1.21) follows similarly, and (1.22) from combining 
(1.20) and (1.21). □ 

It is also convenient to record a version of (1.21) which holds for real powers of 
Ay in the tangential pseudodifferential calculus. Let Fg = (1 + AyY^"^. 



Lemma 1.7. For all s G M, 



[□,n]y_« = P£), + -Q, 

X 

Y_,[n,Y,] = P'D,+ '^Q', 



with P,P' e C°°([0,e);*-i(y)) and Q,Q' G C°°([0, e); ^^(y)) and where Q,Q' 
annihilate constants at x = 0. 

2. Mapping properties of the Laplacian 

By definition, A is symmetric as an operator on C^{X) with respect to the 
Riemannian volume form 

(2.1) dg = x"--'^dxdh{x). 

To keep track of the weighted and Sobolev spaces which necessarily appear here, 
we shall refer all weights to the intrinsic boundary weights. These correspond to 
a non-vanishing positive smooth density on X° with is of "logarithmic" form near 
the boundary 

(2.2) Ub = a—dho, < a e C°°(X) near dX. 
Thus we set 

(2.3) Ll{X) = |« e Ll,{X°y, \u\\b < ooj . 
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Since the conic metric volume form is dg = x^a'uh with < a' e C°°{X), the metric 
L"^ space is 

(2.4) Ll{X) = |« e Ll,{Xy, l^\u\% < ooj = x-^Ll{X). 

From (1.14) it follows that 
(2.5) 



A = x' 



-2 



(a;£),)2 - !t_Jl±IlD, - i{n - 2)xD, + At, 



e x-2Diff^(X). 



Here, T)\E'^{X) is the filtered algebra of differential operators on X which is the 
enveloping algebra of the Lie algebra Vb(^) of all smooth vector fields on X which 
are tangent to the boundary. 

The weighted b-Sobolev spaces x'pH\^{X) are essentially defined by the mapping 
properties of these b-differential, and the corresponding b-pseudodifferential, oper- 
ators. They may also be defined directly using the Mellin transform. Any weighted 
b-differential operator defines a continuous linear map: 

(2.6) P&x'' Difr^(X) =^ P : xPHl{X) — > xP+'' hI'"" {X) V p, Z e R. 

An elliptic b-differential operator, i.e. one for which the characteristic polynomial 
in tangential vector fields is invertible off the zero section of ^T*X, has the inverse 
property with respect to regularity 

(2.7) Pex'' DiS^{X), x-^'P elliptic , u e xPH{;°°{X), Pu e xP+''Hl-"'{X) 

=^ue xPHl{X). 

(See [17] for a detailed discussion of b-differential and -pseudodifferential operators.) 

Such an elliptic operator is Fredholm as an operator (2.6) for all but a discrete 
set of p e R. These correspond to the indicial roots, those values of the complex 
parameter for which the indicial operator is not invertible. The indicial operator is 

defined in general by 

(2.8) P e x'' Diff;"(X) =^ P{x''v) = x'^+^'iliP, s)v + 0{x)). 

This definition depends on the differential, at the boundary, of the defining function 
chosen. Rather than carry the normal bundle information to make this invariant 

we shall simply choose x to be the defining function in (1.3). 

From (1.14), the indicial family of the Laplacian at a boundary component Y is 

(2.9) I{A,s) ^ Ay -i{n-2)s + s'^. 

If = Ao < Ai < A2 < A3 < . . . is the sequence of eigenvalues of Ay, repeated with 
multiplicity, then the indicial roots of the Laplacian are 

(2.10) sf = i^^^±y{n-2r + Xl 
In general, assuming P G x'^DiS^{X) to be elliptic. 



(2.11) P is Fredholm as an operator (2.6) 

—p is not the imaginary part of an indicial root. 
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Thus, and — n + 2 are always singular values of p for the Laplacian. For n > 2 
there is a gap in between: 

(2.12) (— n + 2, 0) is free of singular values. 

Such a gap corresponds to boundary regularity for solutions. Thus 

(2.13) Pex'' DiffJ^(X) elliptic, u e x^HliX), Pu e a;«+''i?^-™(X) 

u e x'^HI{X) provided (p, q] is free of singular values. 

The conclusion of (2.13) (where we assume that p < g to avoid triviality) does 
not follow if {p,q] contains a singular value. For our purposes it is enough to 
consider the special case that 

(2.14) 

p* e {p, q) is the unique singular value in (p, q) and p*<p+l, q < p* + 1. 

In fact let us further suppose that there is only one singular value s € C such that 
p* = — Ims, and I{P,s) is not invertible. Then 

(2.15) u e xPH^X), Pu e x'+ml-"^ ==^u = u' + u", 

Pu' e C°°(X), u" e x'^H^iX), u' = x''{logxyvjX + u' and 

0<j<fe-l 

u' e x''Hl{X) where I{P,xD^) ^ x''{\ogxyvj j = 0. 

\0<j<fe-l / 

Here, k is the order of s as a pole of I{P, s)~^ and x = xi^) is a cutoff which is 
identically 1 near the boundary. If there are several values of s with Im s = — p at 
which I{P,s)~^ is singular then it is only necessary to add corresponding sums to 
(2.15). 

3. Domains and powers 

Applying the general results above for b- differential operators to the Laplacian 

we find 

Proposition 3.1. If n> A then 

(3.1) Dom(Z\) = {w G x'"Ll{X)-Au G ^^(X)} 
is independent of w in the range 

Tl 

(3.2) -n + 2<w<--+2. 

If n = 3 the same is true for w in the range —l<w<0 and for n = 2 

(3.3) Dom{A) = {u e Ll{X);u = c + u', c G C, u'€x'"LI{X), Z\w' G Lg(X)} 

is independent of w for w > sufficiently small. In all cases, A is an unbounded 
self-adjoint operator 

(3.4) A : Doni(Z\) L^iX) 

and Dom(zi) coincides with the domain of the Friedrichs extension; if n> 3, then 
A is essentially self-adjoint. 
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Proof. For n > 3 the constancy of Dom(Z\) in terms of w follows from (2.12) and 
(2.13). Thus, by hypothesis, u E .x"'L^(X), with w in the gap (— n + 2,0) and 
Pu = x^Au e a;~5+^L^. If n > 4 then -f + 2 also lies in the gap, (2.12), so by 
(2.7) and (2.13) 

n > 4, u e Dom(zi) m e a;-^+^i7^(X) 

so then 

(3.5) Dom(/i) = x'^-^Hl{X), n > 4. 

For n = 4, — 2: + 2 = 0is the top of the gap whereas for n = 3 the gap is (—1,0). 
In these cases we deduce only that 

(3.6) Dom(zi) = J u e Pi x'^HliX) : Ziu e i^(X) i , n = 3, 4. 

I tu<0 J 

On the other hand, for n = 2 there is no gap. The hypothesis Au e -^g(^) is 
equivalent to Pu = x^Au e xLg{X). The collapsed gap, at p = 0, corresponds to a 
double root of 7(P, s) = Ay + s^. Thus, (2.15) becomes 

(3.7) u e x-'^H^iX), Au e Ll{X) =^ 

u = c + c' log x + u", u" e x^H^ {X), e > sufficiently small. 

In particular e must be smaller than the smallest non-zero eigenvalue of Ay, 
Xi{Ay). The hypothesis (3.3) on u e Dom(Z\) is therefore just the vanishing of 
the coefBcient c' in (3.7). It follows that (3.3) is also independent of w provided 

W < Xl{Ay). 

We now demonstrate selfadjointness by showing that the unbounded operators 
with the domains described above coincide with the Priedrichs extension of A. By 
definition, A is associated to the Dirichlet form 

(3.8) F{u,v)= I {du,dv)gdg , u,v eC^{X°) 

Jx 

where dg is the metric volume form (2.1). The inner product in (3.8) is that induced, 
by duality, by the metric on T*X°. Following Friedrichs we define 

(3.9) V = Dom(Z\5) = cl{c~(X°) w.r.t. F{u,u) + HuHis} , 

whenever X is a compact conic manifold with boundary of dimension n > 2. Then 
the Friedrichs extension of A is the unbounded operator with domain 

(3.10) Dom(/iFr) = {ueV;Aue L^iX)} . 

where A is the bounded operator "D — > V and L'^{X) C 'D' is a well-defined 
subspace since V C L^^X) is dense. 

The space Dom(Z\2) is independent of which conic metric on X is used to define 
it, since different conic metrics give equivalent norms in (3.8). Moreover it can be 
localized by use of a partition of unity. Thus it is the same, locally, as the same 
space on a compact manifold blown up at a single point. This is well-known and 
easily leads to the characterization 



(3.11) Dom(Z\5) 



\x-^+^Hl{X)+p{x)t> n = 2, 
\x-i+^Hl{X), n>3. 
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where p{x) is a cutoff near the boundary in terms of a boundary defining func- 
tion and V is the analogous domain of the 1-dimensional operator D^, acting 
on a;~^L^([0, oo)). Using Paley's theorem it is straightforward to characterize 
"D C x~^L^([0,oo)) in terms of the Mellin transform: 



(3.12) ueV x-^L(,{[0, oo), 6>0, and 

/dx 
x^'^u{x)— is holomorphic in — oo < 1ms < with 



X 

suMis)\'^dRes uniformly bounded in r e (— oo,0). 

Im s—r 



In the definition of the Friedrichs domain, (3.10), the action of A is distributional. 
Thus from (3.1) and (3.2) it follows that 

(3.13) Dom(ZiFr) = {u€ x-^+'^Hl{X); Au € x-^LI{X)} = Dom(/i), n > 2. 

For n = 2 the argument only needs to be modified slightly. It follows directly 
from (3.3) that Dom(Z\) C Dom(Z\Fr)- Moreover (2.15) shows that u G Dom(Z\Fr) 
has an expansion as in (3.3) except for the possibility of a logarithmic term. This 
however is excluded by (3.12) since it would correspond to a double pole of the 
Mellin transform at s = 0. 

Thus in all case we have shown that the Friedrichs extension has domain Dom(Zi) 
as given in Proposition 3.1. □ 

We also need to describe the domains of the complex powers of A. For integral 
powers it is straightforward to do so. Since 

(3.14) I{A'', s) = I{A, s)I{A, s + 2i)... I{A, s + 2{k- 

it follows that the singular values of A'' are just the unions of the shifts of those 

of Z\, with the possibility of accidental multiplicity to be borne in mind. We are 
particularly interested in the domains of the small powers with real part up to n/4. 
Note that it follows from arguments directly analogous to those above that 

(3.15) Dom(ZiP) = x-^+23'i?2P(x) for p < ^, 2p e N. 

For later applications we need to find the largest real p for which this remains true. 
To do so it is convenient to use complex interpolation. 



Lemma 3.2. For real < p < n/4: the identification (3.15) remains true. The 
domain of z\"/"* is independent of the conic metric defining the Laplacian and is 
given explicitly by 

(3.16) Dom(ZV?) = {X) + p{x)Vn 

where Vn reduces to (3.12) in case n = 2 and in general is defined by 

(3.17) u&Vn <s=^ u e x~^Ll{[{), oo)), 8>Q and 

/dx 
x^^u{x) — is holomorphic in — oo < Ims < with 
x 

I (|s|^ + l)'2~^|suM(s)pc?Res uniformly bounded inr € (— oo,0). 
ill 



Im s=r 
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There exists Sq > 0, depending on the metric and manifold, such that for all 6 € 
{0,So), 

(3.18) Dom(zit+5) = x'^^H^^^\x) + C 
and 

(3.19) Dom(zif) D x-^+'^PH^P{X), V p > n/4; 
by duality there is a restriction map 

(3.20) Doni(zi-P) — > x-^-'^PH-^P{X). 

The map (3.20) is not injective for p > ^ since the inclusion (3.19) does not then 
have dense range, as follows from the computation of the domain of Zi"/^. 

Proof. We first use complex interpolation from the characterization of the domain 
of zA* above. For two Banach spaces X and Y, let [X, FJg denote the complex 
interpolation space at parameter 6 as discussed, for example, in §1.4 of [23]. For 
any positive sclf-adjoint operator the complex powers satisfy complex interpolation 
in the sense that 

(3.21) \}owl{A^) = \DowL{A).Ll\g. 

Furthermore the usual arguments with Sobolev spaces show that the weighted b- 
Sobolev spaces exhibit the same interpolation property 

[x'H^iX), Ll{X)]e = x''H^'{X), < 9 < 1. 

Applying this to (3.15) we conclude that it remains true for p smaller than the 
greatest half- integer smaller than n/4. For n > 2 we may apply the same argument 
again by noting that 

Dom(Z\«) ^ {ue Dom(Z\); An € Dom(Z\«"^)} , 1 < g < n/4. 

This proves (3.15) for all real p < n/4. 

In fact essentially the same method applies to Dom(Z\T ) since we have computed 
Dom(Z\?"i). The condit ion on the McUin transform in (3.17) just represents decay 
at infinity like |s|~"/^ in a uniform sense except for the single factor of s which 
allows more general behavior at s = 0. 

The final characterization of the domains in (3.16) now follows by use of the 
Mellin transform near the boundary and reduces to the same argument as in the 
two-dimensional case. Thus 

(3.22) Dom(ZiT) = |m g x'"'''^ {X), An G i?|"^(X)| 

from which (3.16) follows as before. 

For p > n/4, the fact that Dom(ZiP) D x-^+'^phI^{X) follows from (2.15) and 
(3.18) follows by similar arguments as those above, with the upper bound on 5 
arising from the first positive eigenvalue of the boundary Laplacian. □ 

The identification, (3.18), of the domain of the powers just larger than ^ gives 
a convenient mapping property independent of dimension and of the conic metric 
involved 

(3.23) A ■ x^^H^^^\x) + C ^ x-^+^^H^~^^^\x), S>0 small 
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which has null space C and closed range which is a complement to C. Indeed A — X 
also defines such a map for any A S C and this map is an isomorphism if A is not 
in the spectrum. 

4. Statement of the theorems 

Let £s denote the energy space, of order s e R, of Cauchy data for the wave 
equation: 

£s=Vs(S Vs-i where Vs = Dom(Z\^/^) 

as discussed above. If O is an open set containing a component of dX, we will denote 
by Vg{0) the corresponding local space which is well defined since reduces to 
locally away from the boundary. define £s{0) similarly. 
The Cauchy problem for the wave equation 

(4.1) _ A)u{t) = 0, 

(4.2) w(0) = uo, Dtu{Q) = ui 
has a unique solution 

(4.3) u e C"(M;2?,,) n C^(R;X>5_i) 

for all [uq.ui) G Eg- Similarly the inhomogeneous forcing problem 

(4.4) {Dj - A)v{t) =f,fe C-°^(R;P,_i), / = in t < 

has a unique forward solution v G C~°°(Ili; Pg) with v = in t < 0. 

If the conic manifold is complete rather than compact, away from the conic ends, 
then a similar result holds. In the general case of a manifold with conic ends the 
wave equation (4.1) has a unique solution with compact support in a finite interval 
[— T, T] in place of M in (4.3) provided the initial data has compact support. Since 
the results below arc all local near the boundary and the general case can be reduced 
to this one, for simplicity of presentation we consider only the case of a compact 
conic manifold. 

Definition 4.1. An admissible solution to (4.1) is one of the form (4.3) for some 
s e M with the equation holding in C~°°(M; T>t) for some t gM.. 

We deal here with admissible; solutions to the wave equation, i.e. solutions corre- 
sponding to the Friedrichs realization of the Laplacian. 

The time-translation invariance of the wave equation means that if u is an admis- 
sible solution and e £ C~°°(M) then e{t) * u is also an admissible solution. We shall 
use this below to decompose solutions into positive and negative parts by choosing 
a decomposition 

(4.5) 6{t) = e+{t) + e-{t) + e^{t), e± G S'{R), 

eoc e cS(M), WF(e±) = {(0,±oo)} C S*R. 

The corresponding decomposition of an admissible solution u is then 

(4.6) ii = u+ + M_ + Uoc, Uoc € C°°(]R; I?oc), u±=e±{t)*u, 

where all three terms are admissible solutions. Typically we choose e+ and e_ to 
have Fourier transforms supported in [1, oo) and (— oo, —1]. 
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Global regularity theory for the wave equation shows that the strongly continuous 
group of bounded operators defined by (4.1) - (4.3) 

(4.7) U{t) :Vi®Ll^Vi®Ll 
satisfies 

(4.8) U{t) -Es — > £s V s e ]R, 

by continuous extension for s < 0. The solution to (4.4) is then given by Duhamel's 
principle. 

It is useful at various point in the discussion below to change the degree of 
regularity of an admissible solution; this can always be accomplished by convolution 
in t. 

Definition 4.2. Let ©g be the operator on M given by 

(4.9) K{Qs){t,t')=^{t-t')K{W){t,t') 

where ^{t) is a smooth function of compact support, equal to one near t = Q and 
K denotes Schwartz kernel. 

Lemma 4.3. Let u <E C{M.\£r) he a solution to the wave equation. Then ©gW S 

C{M^;£r-s) is also a solution, and for all s G K, 

QsQ-sU = u mod C'^(R;foo)- 

Away from the boundary the wave operator D^—A is smooth with principal sym- 
bol — \C,W at the point {t, z; r, C) in terms of the canonical coordinates associated 
to the coordinates (t, z). Hormandcr's theorem on the propagation of wavcfront set 
for operators of real principal type therefore applies and shows that microlocal reg- 
ularity, in terms of the wavefront set relative to Sobolev spaces, is contained in the 
characteristic variety and is constant along the null bicharacteristic which foliate it. 
Thus the regularity of a solution at any point of T*(RxX°)\0 is readily describable 
in terms of the regularity of the initial data unless the null bicharacteristic through 
the point hits the boundary at some intervening time. 

We first state a diffractive theorem, which simply says that if there are no in- 
coming singularities at the boundary component Y at some time i then there are 
no outgoing singularities arising at that time. Moreover this regularity can be 
microlocalizcd in r, the dual variable to t. 

Theorem 4.4 (Diffractive regularity). If u is an admissible solution to the conic 
wave equation and for some small e > 0, s G R, i{t,Y) (1 WF^{u) = 0, then 
o(^' ^) <^ WF"(m) = 0. // R}{i, Y) n WF*(u) = 0, then for some open sets l3i 
andODY, u & C{I]V,{0)). 

By time-reversibility of the equation, this theorem implies Theorem L2. It is 
proved in Section 11. 

In order to state the geometric theorem, we need to introduce a second-microlocal 
condition on the incoming singularities at time t = i. 

Definition 4.5 (Nonfocusing conditions). We say that an admissible solution u to 
the conic wave equation satisfies the nonfocusing condition at time i and at 
boundary component Y, with background regularity r and relative regularity I, if 
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for some positive integer N and small e > 

(4.10) WF''(u) n R^jit, F) = and 

(4.11) WF''+' ((1 + Ay)-^u) n R%j{i, Y) = 0. 
If, for some N = N{k) e N, 

(4.12) WF"- {{xD^ + {t- t)Dt)P{l + Ayy'^u) n R%^j{t, Y) = <D, p<k, 

we say that u satisfies the radial regularity condition to order k (and with 
regularity r). The combination of (4.10) with (4.12), but for regularity r + I: 

(4.13) WF''+' {{xD^ + {t~ t)DtY{l + Avr^u) n i?ij(t, F) = 0, p < fc, 

will be called the conormal nonfocusing condition (to order k with background 
regularity r and relative regularity /). 

These conditions measure the extent to which the solution contains a wave col- 
lapsing radially onto the boundary, up to relative regularity order I. The strongest, 
conormal, version asserts that after smoothing u in the tangential variables a conor- 
mal estimate at the surface x = t — t holds to order k. As an important example, 
the fundamental solution s\iit\fA/ \/~A with pole close to the boundary satisfies the 
conormal nonfocusing condition with r < — n/2 + 1 and I < (n — l)/2 for any p (see 
§14). 

Theorem 4.6 (Geometric propagation). Let u be an admissible solution to the 
conic wave equation. If p € R^_j_Q{t,Y), r'^(p) n WF'' (w) = for small e > 
and u satisfi,es the nonfocusing condition of Definition 4.5 with regularity r + 1 > r', 
thenp(^ WF' '(m). 

This theorem represents a sharpening of Theorem 1.3. We prove a weaker version 
of it in Section 15 and then obtain the full theorem in Section 16. 

The "edge structure" on the product x X is discussed in the next section 
and in particular the related scale of weighted edge Sobolev spaces is defined there. 
These spaces are appropriate for the description of the boundary regularity of ad- 
missible solutions to the wave equation. A crucial role in our proof of the geometric 
propagation theorem above is played by the following result on decay relative to 
weighted edge Sobolev spaces. 

Theorem 4.7 (Division theorem). If u is an admissible solution to the conic wave 
equation satisfying the nonfocusing condition in Definition 4-5 with r + I < n/2 
then there are open sets I 3 t and O dY in X such that 

(4.14) u e a;''+'-"/^ii;+'-'=(/ x O). 

If u satisfies the conormal nonfocusing condition to order k = 1 and if r + 1 < 1 
then in addition 

Dtu e a;'-+'-"/2-l/2-e^r+i-ft-l^j X O) V e > 0. 
We prove this theorem in Section 14. 

One can use energy conservation and the diffractive theorem to show that any 
solution to (4.1)-(4.3) is in a;*~"/^i7^ as long as —n/2 < s < n/2 and that any 
solution satisfying WF''(u) n i?±.j(f) = is in .x''-"/2i/J(]R x X), locally in time 
near t = t, &s long as —n/2 < r < n/2. The nonfocusing condition thus leads to a 
stronger growth estimate than is given by energy estimates alone. 
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The fact that the division theorem yields a stronger result in the presence of a 
(first-order) conormality assumption is closely related to the fact that if a solution 
to the wave equation is conormal with respect to the incoming surface x — t = x in 
t < it is conormal with respect to the corresponding outgoing surface x = x + t in 
t > 0. This follows in turn from the corresponding result on the Cauchy problem. 

Theorem 4.8 (Conservation of conormality). Let (uq^ui) (E £s be conormal with 
respect to the hypersurface {x = x} for some x > sufficiently small, then inx > 0, 
the solution to (4.1)-(4.2) is conormal with respect to 

{x = \x-t\}U{x = x + t} 

forO <t < 2x. 

The proof of this theorem appears in Section 13. 

A microlocalized version of Theorem 4.8 allows us to show the conormality of 
the diffracted front, subject to radial regularity of the incident wave. 

Theorem 4.9 (Conormality of the diffracted front). Let u be an admissible solu- 
tion to the conic wave equation satisfying the radial regularity condition in Defi- 
nition 4-5 at i,Y to every order k & Nq and suppose that for some small e > 
and some p S RQ{t,Y), T'^{p) (1 WF{u) = then, microlocally near the outgoing 
bicharacteristic through p, u is conormal with respect to x = t — i. 

The proof of this theorem is in Section 15. As a corollary of this result we may 
refine Theorem 4.6, concluding that if in addition the full conormal nonfocusing 
condition, (4.13), holds for all k, then the solution is conormal to the surface x = t—t 
near p. 

Although it follows from the results above, we nonetheless restate part of the 
result on the fundamental solution discussed in the introduction; in fact, the proof 
of this result, which occurs Section 16, is crucial in the proof of the full version of 
Theorem 4.6. 

Theorem 4.10 (Fundamental solution). Let Em be the fundamental solution to 
the conic wave equation with pole atm = {x,y) G X°. If x is sufficiently small then 
Efn is conormal with respect to {x -\- x — t} away from the wave cone emanating 
from fh, and is of Sobolev order ^ — S there for any S > 0. 

This is equivalent to Theorem I.l of the Introduction. 

5. Edge pseudodifferential calculus 

The edge calculus of pseudodifferential operators was introduced by Mazzco [14] 
as a class of operators on any compact manifold M with boundary having a fibration 
dM, TT : DM N with fiber F. In this paper, M = M x X, A'' = R, and tt is the 
product fibration 

TT : R X 5X ^ E, 7r(i,p) = t. 
The noncompactness of M in the situation at hand necessitates only minor changes 
to the calculus, namely keeping supports proper. Our operators generally have 
compactly supported kernels. In this section, we discuss the edge calculus in the 
general setting, as this involves no increase in complication over the special case of 
R X X. Although there is no treatment in [14] of edge microsupport or wavefront 
set, the properties of these objects follow easily from the properties of the calculus 
discussed in [14] much as the properties of the conventional wavefront set and 



22 



RICHARD MELROSE AND JARED WUNSCH 



microsupport follow from the properties of the pseudodifferential calculus on closed 
manifolds. First wc give a brief synopsis of the edge calculus and its properties. 
A Lie algebra of C°° vector fields Vc(M), associated to tt, is given by 

V GVe <^=^ V is tangent to the fibers of n at dM. 

If {x, y, z) are local coordinates with x a defining function for dM, y coordinates 
on TV lifted and extended to functions on M, and z restricting to coordinates on 
the fibers, then Ve(M) is locally spanned over C°°(M) by 

xDx, xDy and D^. 

Thus, there exists a vector bundle '^TM (the edge tangent bundle) such that 

Ve(M) =C°°{M]^TM). 

Let «'T*M (the edge cotangent bundle) denote the dual of ^TM; sections of ''T*M 
are then locally spanned over C°°{M) by dx/x, dy/x and dz. By the edge cosphere 
bundle we mean the quotient 

''Snd = (°T*M\0) /M+. 

There is a canonical bundle map °T*M T*M since edge vector fields are 
smooth up to dM. This map is an isomorphism over M° so over the interior °T*M is 
a symplectic manifold, with symplectic form given as usual by the exterior derivative 
of its canonical one-form; this form becomes singular at the boundary. 

For fc G N, let DiffJ(M) be the space of differential operators spanned over 
C°°{M) by operators Vi . . .Vm, Vi E Vc(M), m < k. There exists a canonical 
(principal) symbol map, "a^, associating to P € Diff™(M) the polynomial function 
on the fibers of ^T*M extending the usual symbol map over the interior; it gives a 
surjective map 

''(7m ■ Diff™(M) — > homogeneous polynomials of degree m on ''T*M 

with null space precisely Diff™~^(M). 

In the particular case of interest for this paper, Ve(K x X) is locally spanned by 
the vector fields xD^, xDt, Dy over C°°(R x X). If we let 

^dx ^ dt 
^—+X—+Vdy 

x X 

be the canonical one- form on ''T*(M x X), then naturally 

Vi(a;£>^) = ^, Vi(xA) = A, = m 

are dual coordinates on the fibers of ''T*M. 

The edge calculus of pseudodifferential operators, defined in [14], arises as a 
microlocalization of Diff*(M). Let C°°{M) denote the space of smooth functions 
on M vanishing, with all derivatives, to infinite order at dM, and let C~°°(M) be 
the dual to the corresponding space of densities, C°°(M; Q). The space \E'J(M) is a 
graded algebra of operators on C~°°{M); as we will frequently use weighted edge 
operators, we will depart from the notation of [14] in carrying along the weight as 
an index in the calculus. Thus the bigraded space of operators ^^'\M) enjoys the 
following properties: 

• ^f™'' (M) is a graded *-algebra. 

• ^'^''(M) = a;'*™'°(M) (and the latter space is the space denoted *™(M) in 
[14]). 

• x^ Diff™(M) c *™''(M) for all m e N, / e Z. 
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• The maps '^am extend to 

V™j : *r''(M) x' [s;i^{^T*M)/S;i-^'CT*M) 

the range space for ®a can be conveniently identified with C°°{^S*M). 

• The symbol map is a homomorphism of *-algebras. 

• The sequence 

— > *™-i''(M) — > ^^'\M) — > [s'^g(<'T*M)/S'™-i(''T*M) 

is exact and multiplicative. 

• If A e *™''(M) and B e *^''''(M) then 

t 

where the Poisson bracket is computed with respect to the singular symplectic 
structure on '^T*M described above. 

• ''a(x-^Ax^) = °(j(A) for all A e ^^(M), ^ € M. 

In the case M = R x X, the elements of \1/™''(M x X) may be represented locally 
in the form 

(5.1) Au{t, X, y) = (27r)-"-2a;' j e'(^-i)f+(J'-s'')''+*^^ 

h{t, X, y, s, y', T, ^, r], X)u{xs, y', t - xT)x^+''ds dy' dT, 

where 6 is a classical (polyhomogcncous) symbol of order m with ^, r], A as fiber 
variables; '^amiA) is then the equivalence class of b. 

An operator A € *^''(M) is said to be elliptic at p e *=S'*M if a{A) has an 
inverse in x-^[S~^{^T*M)/S-^-\''T*M)], locally near p. Assuming that the or- 
ders are clear we may suppress indices and so p G Elle(A). There is a related notion 
of microsupport for edge pseudodifferential operators, corresponding to the fact 
that the composition of operators gives an asymptotically local formula for the 
amplitude, b in (5.1). If A e *™''(M), the microsupport of A, WF^(A), is the 
closed subset of^S*M given locally by ess supp b, the conic support of 6; it has the 
following properties: 

. wf;(ab) c WF'^{A) n WFi(B). 

• WF'^lx-'^Ax'') = WF'^{A) for all s e M. 

• If A e *™''(M) and p G Elle(A) there exists Q e "'(M) such that 

p i ^F'^^QA - J) U WF;(AQ - J). 

• If A e *™''(M) and WF;(A) = then A e ^'-°°''(M)— note that this is not 
a totally residual operator. 

There is a continuous quantization map (by no means unique) 
Ope : a;'5^g(^T*M) ^ *^''(M) 

which satisfies 

V™,KOpe(a)) = [a] e a;'5™hg(^T*M)/S™-i(^T*M) V a e x' S^^^i^T* M) and 

WFgOpe(a) C ess supp(a). 
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Associated with the edge calculus there is a scale of Sobolev spaces. For integral 
order these may be defined directly. Thus for A; e N and any s e R we set 

(5.2) H^'^iR xX) = {ue x^Lli,,{R x X)- 

{xDtYPk-rU e x'Lli^^{R X X) V Pk-r e Difr^-''(X) and < r < A;}, k€ N. 
For negative integral orders we can similarly define 

(5.3) for k e -N, H^''{R x X) 3 u ^ u G C-°°(R x X) 

and 3 Ui,r e x'Lli^^{R x X), Pk-r,i e Difr^-''(X), z = 1, . . . , TV, 

N k 

with W = ^ y^^{xDtY Pk-r,iUi,r- 
i=l r=0 

For general orders, the edge Sobolev spaces can be defined using the calculus. 
Definition 5.1. u e H]r\M) ^ *™'-'(M) • u c Ll{M). 

Note that we have chosen to weight these Sobolev spaces with respect to the b- 
weight, not the metric weight. Note also the change of sign on I. Since H^'\M) = 
x^Hf^^iX), wc will often use the notation H^{M) = H;r-"{M) and write the 
x-weight explicitly. In the case of interest in this paper, when M = M. x X is 
noncompact, we will consider only edge Sobolev spaces local in t, without writing 
this explicitly. The corresponding i^-based edge wavefront set plays a fundamental 
result below. 

Definition 5.2. If m e H-°°'\M) then ^¥'^'\u) C "S*!^ is defined by the con- 
dition that p ^ WFr''(u) iff there exists A € *™ "'(M), elliptic at p, such that 
Au e Ll{M). 

The usual properties carry over to these spaces: 

• WFr''(M) is closed. 

• For all I e M, W¥'T'\u) n ^=5*7^° = WF™(u) (recall that «S'*M° and S*M° 
are canonically isomorphic). 

• H^'\M) = C°°{M), [jra,i Hr\M) = C— (M). 

• Complex interpolation holds: 

• If A e *™''(M) then A : Hf'^'{M) — > if™'-™'''+'(M). 

• For m < m', 

wFr''(u) c wFr'''(w). 

• If w G H-°°'''{M) and A G *^''(M) then 

WFT'^+^' {Au) c wf;(A) n WFr+'='''(w). 

• If u e ff-°°'''(M) and A e *^''(M) then 

WFr+'='''(u)\WFr''+''(Au) C (EUe^)"^. 
We now prove a less standard result. 
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Proposition 5.3. //we i/~°°''(M) and I' < I then for all 6 e [0, 1] and m, m' e 

K, 

^pem+(i-e)m',ei+(i-e)j'(^^ ^ WF™''''(u) n WFr''(w). 

Proo/. Suppose p ^ WF™''''(u) n WF™''(m). Then there exists A! G 4-™' "''(M) 
and A e *™'~'(M) both elhptic at such that Au, A'u e ^""o™ elHptic 
regularity it follows that if i? G 'J'"'"(M) has WF^ concentrated near p then Bu G 
if™''(M) DH^'- (M). By interpolation of weighted edge Sobolev spaces it follows 
that Bu G TjfC^-^O.^e-'Oj-^^) j^j^^ ^i^g j.ggyi^ follows, □ 

We will require some results about edge-regularity of solutions to the wave equa- 
tion. 

Proposition 5.4. For all p G 1R+, 

x^H^^iR X X) = Hl^iW; LI{X)) n x^Ll^{R- Hl{X)) and 
x-^H-^{R X X) = H-/^{R; LI{X)) + x-''LUR; H-^{X)). 

Proof. All spaces are local in t (by definition); multiplying by any (/> G C^(M) it 
suffices to assume that supports are compact. Fourier transformation in t coupled 
with an interpolation argument shows that for p G N 

HlSR■Ll{X))r^xPLlM■,Hl{X)) = n ^'K'A^-^HiiX)). 

The latter space is equal to xPHp{R x X) as defined in (5.2), proving the first result 
for p G N; the second follows by duality. The results for general p G R follow by 
interpolation. □ 

Proposition 5.5. Ifu is a solution to the wave equation in L'^g^{R;'Ds), with \s\ < 
n/2, then u G x"if|""^^(M x X). 

Proof. Certainly u G H^^^{M.;T>s-2r) for all r G Z. So, by interpolation, this holds 
for all r G M. By Lemma 3.2, u G Lf^^{M.;x^-''/'^H^{X)) n iJf„^(R; Lg(X)), hence 
the result follows from Proposition 5.4. □ 

6. BiCHARACTERISTIC FLOW 

The canonical one-form on '*T*(M x X) is 

, dt ^dx 
\—+C—+Vdy 

X X 

hence the symplectic form is 

dX Adt d£ A dx Xdx Adt 
u) = 1 ^ \-dr]A dy. 

X x x'^ 

We can now write the symbol of the d'Alembertian 

p = a{n) = X^/x^ - g{x, y, ^dx/x + T]-dy) = 

and 

Po = cr[V^ - Aq) = 2 
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where Aq is defined by (1.15). Let E = {p = 0} C ''T*{R x X) denote the 
characteristic variety of the d'Alcmbcrtian. 

Let Hg and Hg^^ denote the respective Hamilton vector fields of p and po on 
''T*(IR X X), near a boundary component Y of X. Thus 

2 2 

(6.1) ^Hg = ^Hg, +W = HY + {e + ho{y, v))di + X^dx + ^xd, - Xxdt + W, 

where W is the Hamilton vector field of p — po, hence 
(6.2) 

and where i/y is the Hamilton vector field in {y,r]) for {l/2)ho{y,r]), i.e. is the 
geodesic spray in Y. 

Note that Hg, Hg„, and W are all homogeneous of degree 1 in °T*(R x X), and 
that (x^/2)VF is less singular than (x^/2)iJg at a; = 0, and vanishes at ry = 0. Thus 
if T*(R X X) denotes the fiberwise radial compactification of °r*(M x X), 

(6.3) {xy2)Hg, {xy2)Hg„ eh{7j)iVbCT* (R x X)\0), 

(6.4) {x'^/2)W e .T/i(r/)^Vb(°r*(R x X)\0). 

Note also that the vector field Hg is tangent to the incoming and outgoing sets 
R±.i{Y), i?±.o(y)) which we now regard (by homogeneity) as subsets of ®S'*(R x 
X). These incoming and outgoing manifolds arc the interiors of smooth manifolds 
with boundary in '^S'*(]R x X), and we define their boundaries as follows, with IC 
and OG standing for incoming and outgoing manifolds respectively. 

Definition 6.1. Let 



\C±{i,y) = Rl^j{i,y) n ^5^x9x(IR x X) 

OG±{t,y) = i?|_o(f,y) n ^^^^^^(M x X) 

with the same convention for omitted indices as was used for ,(f, y). 

For Y a boundary component of X, we define a map 

T :^T*{Rx X)\Od U ^Y 

constant in the fibers, which is approximately invariant under the flow of Hg and 
which will serve as a useful localizer. The subset U on which T is defined is a conic 
neighborhood of IC. The map T is constructed as follows: consider the data (a;, H) 
of Theorem 1.2 as identifying a neighborhood of y in X with a neighborhood of 
a; = in the model cone X = R_|_ x Y, which we now equip with the model metric 
go. Then for any point q near IC in ^T*{R x X), set 

(6.5) '^{q)= 1™ ^Yexp {s{x'^/2)Hgg), 

where 

and where Try is the projection onto the factor Y of X. As will be shown below, 
the signs are chosen so that T(q) is the limit of the projection on Y of the unique 
geodesic through q, as it heads toward the large end of the model cone X, i.e. the 



Soo = ho{r]{q)) ^ ^(sgn0)^ — arctan^^ , with 9 = 



WAVE EQUATION ON CONIC MANIFOLDS 



27 



point at infinity from which the geodesic emanated. To see that T is well-defined 
and smooth, note that under the flow along {x^ /2)Hg^, ^" = 2^^', hence 

$(s) = Ctan(Cs + 6*) A(.s) =Dscc{Cs + 9) 

(6.6) x{s) = Esec{Cs + 9) t{s) = -Etein{Cs + 0) + F 

hoiv) = G. 

Since ^' = _|. hQ{ri), we compute C = ho{ri{q))i and 6 = a.vcta,n^{q)/ho{r]{q))^ , 
so that 

depends smoothly on q if U is chosen small enough (sgn^^ is constant on compo- 
nents of U), and as s — > Soo, x +qo and is strictly increasing on s € [0,oo); 
simultaneously, t ±oo. Thus, since {y,r]) are undergoing geodesic flow, 

(where Hy is geodesic flow on Y with metric /iq) is manifestly smooth in q. By the 
definition as a limit along flow- lines of Hg^^ , we also have 

Now we turn to the perturbed flow Hg-. Let V denote the rescaled vector field 
(xV2)(A2 + ^2 ^ h{r]))-^/'^Hg on ^S'*(M x X), and let the flow along V be paramet- 
rized by s. Note that under the bicharacteristic flow of {x^ /2)Hg on T*(R x X), 
as \t\ — > oo, s — > C~-^(±7r/2 — 6), hence ^/A — > ±C/D; moreover, |^|, |A| — > oo in 
this limit, while r] remains bounded. Thus we have established: 

Lemma 6.2. Every maximally extended integral curve of V over dX contains in 
its closure exactly one point in IC and one in OG; the former lies over the point 
T(p) e Y , for any p along the integral curve. 

7. Construction of symbols of test operators 

We now write down the symbols of the operators to be used in the commutator 
estimates in §8. As usual, we work in a product neighborhood of a boundary 
component Y as described in Theorem 1.2. For points yi, j/2 S Y, we let d{yi,y2) 
denote distance with respect to the metric ha- 

Let X G C°°(IR) vanish for x < 0, be equal to 1 for x > 1, and be nondecreasing, 
with smooth square root, such that x' also has smooth square root. Choose tp{x) € 
C(i"(K) to be equal to 1 at x = 0, be supported in (—1, 1), with derivative supported 
in (—1, —1/2) U (1/2, 1), and to be the square of a smooth function; let (sgna;)V''(a;) 
also be the square of a smooth function. For positive constants Ci let tpi{x) = 
tp{x/ei) and Xi{x) = xi^-ff-i)- 

First we consider test symbols at incoming radial points. 

Given m,l G R, y G Y, and x G R+, deflne a nonnegative symbol on ^T*X by 
setting 

{<,i,±? = x(±A)x(±Oxi(a:^ - a^o + vt)xi{-x + x^- v't) 
•V2((d(T,y)2-fe)+)V3(M»7)V|A|)X4(t + e4)V'5(M*,ar,t/,A,^,r/)/A2)(±A)™a;'. 
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where v < 1 < v' , xo < x < xi, 5 > 0, and we have written T = T{t,x,y,X,^,r]) 
for the map defined by (6.5). We will assume that xi is sufficiently small that the 
perturbation term W — Hg — Hg„ has the property 

{xy2)W = ho{r])^0{x)d^ + 0{x)d^ + 0{x)dy) 
(from (6.2)) where 

(7.1) all 0{x) terms are bounded by 10~^ when x < xi 
and furthermore that 

(7.2) (a;/2)^^ < 10-2 fora;<a;i. 

ox 

Since (x'^ /2)Hg^d(T ,y)'^ vanishes, (6.4) yields 

(7.3) I {x'^/2)HgdiT, y f\< Axhij]) 5 when x < xi 

for some constant A. 

Observe that on suppa™ ; _|_, both 

\>?-e-h{r,)\ 



A2 

and 

Hence 



(7.4) ^l-e5-e^<|ei/|A|<yrT^ 

so that suppa™ ; _|_ is localized arbitrarily near R±j. 

We choose 5 small enough that 5x\ < £2/2, hence suppV'2(-) fl {x = xi} ^ 0. 
We now choose the other sufficiently small and v,v' sufficiently close to 1, so 
that aj^ ;_-|- has support in an arbitrarily small neighborhood of the closure of a 
bicharacteristic segment 

{t = s, X = X — s, y = y, A = ±1, ^ = ±1, r] = 0;s G [0, x]} 

passing through the point q with coordinates x = x, y = ij, t = 0, ^ = ±1, t] = 
and hitting IC at time x. Note that we may translate the t variable freely without 
changing any of the properties of a™. 
We now evaluate, term by term, 



2 



^s('^m,i,±)^- 



The term containing 

YHgXi{x -xo + vt) = xi(a^ + vt){i -v\ + O{x)ho{ri) 2 )x 

has sign ± on suppa^ ; _|_, since (7.4) imphes |^| > |A|(1 — 65 — €3), hence by (7.1), 
sgn(^ — vX) = sgn^ on suppa'" provided w < 1 — £5 — £3 — 10^^£3. 

Similarly, the term involving [x^ /2)HgXi{—x + Xi — v't) can also be made to 
have a positive derivative if v' is sufficiently greater than 1. Shrinking £3 and £5 as 
necessary, these conditions can always be achieved. 

The term containing ^Hgip2{{d{T,y)'^ — 5x)+) gives 

i;'2{{d{r,yf-Sx)+)-{-5^x + 0{x)h{rj)i) 
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where, by (7.3), the 0{x) stands for something bounded by Ax. Hence if 

(7.5) 5\^\>Ah{rj)K 

this term has sign ±; by (7.4) and the support property of tp3{h{r])i /\X\), the 
condition 

^>^^^ 

V 1 - - 4 

suffices to ensure (7.5). This can be achieved by further shrinking £3 as necessary. 

The term involving {l/2)x'^Hgip3{h{'r])^ /\X\) can be evaluated using (6.1)-(6.2), 
which show that Hgh{r]) = ^xdh/dx and hence 

so by (7.2) this term has sign ±. 

The term involving {l/2)x^HgXi{t + £4), unlike those discussed previously, has 
sign ^. This term is supported in a region in which we will assume microlocal 
regularity: its support is in 

< e4, X G (xo, xi), (i(T, yf < 6x1 + £2, h{rifl^/\X\ < £3}; 

this lies away from dX but inside an arbitrarily small neighborhood in ''S'*(IR x X) 
of an arbitrarily specified point in j. Call this term e, for "error." 

The term involving ^Hgip^ is supported in p/A^ > £5/2, hence vanishes iden- 
tically on the characteristic variety E. Call this term k. The terms arising from 
{l/2)x'^Hgx{±X) and {l/2)x'^ Hgx{±£,) are supported in a compact subset of T*(Mx 
X); let c be their sum. 

Finally, the factor (±A)'".t' has derivative (m+?)^(±A)'", which has sign ± sgnm, 
i.e. ± as long as m + ^ > 0; this term, of course, has the same support as j _(_ 
itself; denote it ±(a')^. 

All the nonncgativc rcsp. nonpositivc terms described above can be arranged 
to be squares of smooth functions resp. minus squares of smooth functions. We 
organize the information gleaned above as follows. Let 

q = {i, X, y, X = ±1, ^ = ±1, r? = 0) 

be a given point in R±j with x sufficiently small; let fl be the closure of the 
bicharacteristic connecting q to the boundary: 

n = {(t + s, x-s, y, A = ±1, ^ = ±1, r) = 0) : s e [0,x]}. 

Thus we have shown: 

Lemma 7.1. Provided m+l > there exists a symbol aj^ ; ^ of order m and weight 
I in''T*{Rx X) such that 

(7.6) Y^,i<i,±f = Hd'f ±Y^b'j+e + c + k 

j 

where a! = a ■ (±(to + OO^i ^^PP'^m i ± arbitrarily small neighborhood offl, 

suppe is an arbitrarily small neighborhood of q, suppc is compact in ®T*(R x X) 
and S n supp(fc) = 0. 
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Next we consider test symbols at outgoing radial points. Let V' be as above. 
Given x, y, set 

(CU)' = x(±A)x(TOxi(-a; + X1+ vt)xi{x - xo - v't) 

■ M{d{r, yf + fe)+)^3(/i(r?) V|A|)X4(e4 - t)V5(p(i, x, y, X, r?)/A2)(±A)™a;' 

where v < 1 < v' and xo < x < xi. As with a™, we can choose constants 
a small enough that supp a"^] ^ lies in a small neighborhood of the closure of a 
bicharacteristic segment 

n = {(t = -s, X = x~ s, y = y, A = ±1, ^ = =f1, V = 0); s e [0, x]} 

passing through the point q with coordinates x = x, y = y, t — 0, ^ = ±1, rj = 
and emanating from d{M. x X) at time —x. Moreover if m + / < 0, all terms in 
{x"^ /2)Hga'^\ _l_ can be arranged to be ± squares of smooth functions, with the 
exception of a compactly supported term, a term supported away from S, and, 
most importantly, the term involving 

^if,(^3(M^)V|A|)). 

Let e denote this error term. Then 

supp(e) C supp(a°;'X±)n {/^(77)V|A| G [63/2,63]}. 

This is a subset of the complement of i?±,o inside any given positive conic neigh- 
borhood of Q.. 

More generally, let 

q = {i, a; = 0, y, A = ±1, ^ = Tl, = 0) 

be any point in OG; let f2 be the closure of the short bicharacteristic extending 
from q to {t = t + x, x = x, y = y, A = ±1, ^ = =f1, rj = 0), so 

n = {{t = t + s, x = s, y, A = ±1, C = ±1, = 0) : s G [0, x]}. 

Lemma 7.2. Provided m + l < there is a symbol _(_ of order m in ''T*(M x X) 
such that 

(7.7) ^H,{aZ,±f = ±{a f ± ^ 6| + e + c + A: 

3 

where supp(a^* _|_) is an arbitrarily small neighborhood of il, a' = a{+{m + l)£,)^ , 
supp(e) is contained in the complement of R±^o in o,n arbitrarily small neighborhood 
of VI, supp(c) is compact and S fl supp(fc) = 0. 

8. Propagation of edge wavefront set 

In this section, we prove a theorem on propagation of singularities for the edge 
wavefront set which is a central ingredient in the diffractive theorem (Theorem 4.4), 
the geometric propagation theorem (Theorem 4.6) and the proof of the conormal 
regularity of the diffracted front in Section 15. 

Theorem 8.1. For u € H~°°'''{I x [0,6)^; x Y), a distributional solution to the 
wave equation \3u = 0, with t € I CM. open, the following four propagation results 
hold. 
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i. Ifpe \C{t,Y), m > 1 + {n- l)/2 and WF'"(u) n R^^ j{t,y{p)) = then 
p^WFT'^'iu) for all I' < I. 

ii. For any m e R, WF^'^u) n ''S^.^Q^iM. x X)\{\C{t) U OG{i)) is a union of 
maximally extended integral curves ofV = {x^ /2)[\^ + ^'^ + ho{rij)~^^^Hg. 

iii. IfUc ""S^^Qxi^ X X) is a neighborhood of p G OG{i,Y) and WFT'^{u) n 
[7 c OG then p ^ WFf ■'(Zif.u) for k G No provided M < m - 2k and 
M <l + {n~l)/2. 

iv. IfR^u e H~°°^^ {R X X), with R given by (1.17), and p & OG{t) has a neigh- 
borhood U c''S;.^y(^x^) such that WFT'\R^'u) nU C OG, for < f < j, 
then p ^ WFf (R^u) for j G Nq provided M < m- j and M < 1 + {n-l)/2. 

Remark 8.2. This theorem correspond to propagation into, within and out of the 
boundary. The first two parts, together with the fact that WF™''('u) is closed, 
can be combined with Lemma 6.2 to conchidc that for p G '^S^^g-^{M. x X)\OG, if 

WF'"(«) n R'j^j{t{p),Tip)) = then p i WF™'''(u), for < I. 

The conclusion of (iii) implies, by closedness of the edge wavefront set and 
Hormander's Theorem, that in fact all of R^ Q{i,y{p)) is absent from WF"^{u). 
This third part of the theorem is trivial, however, when it is applied with I = s — n/2 
to a solution in Vg. Therefore this part of the theorem (and the fourth part likewise) 
is useless in the absence of a "division theorem" yielding better x decay of u than 
is given by energy estimates. 

Proof. First consider (i). By assumption, u G H^'\M. x X), locally near {t} x Y for 
some q gM. and p G \C{t,Y). We shall prove the following statement: 

If m' >l' + n/2 - 1, u e i?-°='''(/ X [0, e) xY), p(^ WFf'^'{u), 

^^'^^ and R±jit,y{p)) nWF'^'+^/^iu) = then p ^ WFf+^^^^''{u). 

To do so, choose 

(8.2) As G *™'.-''-/2+l(K X X) = Ope [M\^,V)aZ',-,-n/2+l,±\, 

where a™ is constructed in Lemma 7.1 and 

^i(A,^,r?)=V((A'+^' + M??))^) 
with ^{x) smooth, equal to 1 for a; < 1/2 and for a; > 1. Choosing the supports 

sufficiently small, we see that the error term e in (7.6) will have support in the 

complement of WF™'+^/2(u). Note that Hgips is supported in {1/(2(5) < +^2 

Kv) < 1/^}- 

Thus if A'g has symbol ipgd' with a' as in (7.6), 

(8.3) [□, A*sAs] = ±{A'sy{A's) ± ^ BljBsj + Es + Ks + Rs + Ss 

j 

where 

• WF;(/r5)nE = 

• A'gG 5'r''^^'"''"^(M X X) with A'g ^ A' = A'^ in ^f^^^''~'''~^{R x X) for 
all e > as (5 ^ 0. 

• EgG -2i'-n(]gx -g bounded in 6 with WF'^{Es) uniformly bounded 
away from dX and contained in (WF^'+^Z^ u)^, 
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• Rs bounded in ^f2™' -2r-n(R ^ X) 

• Ss is bounded in ^f^™ +1.-2; -n^jj ^ X) and, as 5 ^ 0, converges to in 
q,2m'+i+e,-2i'-n(^^ ^ j^^ £qj. all e > (this is the term whose symbol involves 

Applying (8.3) to u and pairing with u with respect to the inner product on L"^ 
yields 

(8.4) WA'^uW] - {Ssu,u)g < \{Esu,u)g\ + \{Ksu,u)g\ + \{R5U,u] 

the integration by parts is justified since u G H^°°''- (M x X) and nA*gAsu G 
jjoo,-l -"■(M.xX). All terms on the right-hand side are bounded uniformly as 5 — > 0. 
A weak convergence argument now shows that |l^ou||^ < oo, hence WF™ ^^-^ (y) = 
is disjoint from the elliptic set of A'q-, the shift by n/ 2 in the x weight here comes 
from the difference between and L^. This proves (8.1). 

Ifg > I + {n — 2) /2, then iterative application of (8.1) proves that p ^ WFT'\u) 
directly. If q < / + (n. — 2)/2, however, a further argument is needed. 

Supposing Iq = sup{r;p ^ WF^'^{u)} < I, we wish to arrive at a contradiction. 
We will employ an interpolation argument illustrated in Figure 1 . We have already 
shown that Iq > q — (n — 2)/2 and by hypothesis, u e i?|''(IR x X) (at least locally), 
so Proposition 5.3 shows that 

(8.5) p ^ y^peg+{i-e)m,ei+{i-e)io (u) y 9 e [0, 1] . 

In particular, since / < m — (n — l)/2 (by the hypothesis of the theorem) and 
g < / + n/2 - 1, 

e' = {m-lo- n/2 + l)/((m - la - n/2 + 1) + {I - q + n/2 - 1)) G [0, 1]. 

If m' = 9'q + (1 - 9')m and = 9'1 + (1 - 9')lo, then m' = I' + n/2 - 1, and by 

(8.5) , p ^ WF"''''~'(u) for all e > 0. Hence applying (8.1) iteratively shows that 
p ^ WF™'' ~^(u), with I' — e > lo, which is the desired contradiction. Thus (i) is 
proved. 

To prove (ii), i.e. to show that regularity propagates across *S'g^(M x X) up to 
(but not including) points in OG, we appeal to the standard proof of Hormandcr's 
theorem on propagation of singularities for operators of real principal type by use 
of positive commutator estimates. This applies microlocally near all characteristic 
points where the rescaled Hamilton vector field V = {x^/2){X^ + + \rif)~^^'^Hg 
is non-zero, hence away from IC U OG. See [16] for an analogous discussion in the 
context of the scattering calculus; note that in that case the propagation result may 
be reduced to Hormandcr's theorem whereas in this case it is an analogue of it. 

Now consider (iii), first for A; = 0. In this case we may simply suppose that 
M = m < I + {n — 1) /2. We proceed much as in the proof of (i) above. The result 
follows by iterative application of the following assertion: 

(8.6) Ifm' <l + n/2-l, ue H-°^^\l x X) and p ^ WFr'''(w) then 



P i '=^Rxax(K X ^) n WFe™'+i/2''(u) \ OG ^ p ^ WFr'+'/'''(«). 

Note that the final hypothesis on p here is equivalent to the existence of a neigh- 
borhood U oip such that U n WFr'+^/^''(u) C OG. 
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Figure 1 . The interpolation argument in part one of Theorem 8.1: 
we begin with global regularity of order {q, I) and microlocal reg- 
ularity of order (m, ?o)- Interpolation gives microlocal regularity 
of order {m',l') and iterative application of (8.1) is used to move 
along the horizontal line and obtain microlocal regularity of order 
(m, I' — e) with I' — e > lo- 



To prove (8.6), choose e ^f' ' "/^+^(MxX) as in (8.2) with ajj^, _;_„/2+i,± 
replaced by a^*_;_„/2+i ±' where a°"* is constructed in Lemma 7.2 and supports 
are kept small, corresponding to the implicit neighborhood, U, in (8.6) 
symbol ips<i' with a' as in (7.7), then 

[□, A*sAs] = ±{A'sy{A's) ± Bh + Es + Ks + Rg + Ss 



If A'g has 



where 



WF'^iKs) n E = 

- rn' + 'i-, — 2/— n 



A', e *e 



A'q in * 



X X) for 



X X) and A'g A' 
all e > as (5 0. 

Eg e 4r2m'+i,-2i-n(]g ^ x) and WF^iEg) C U\OGit), uniformly in S. 
Rg is uniformly bounded in ^/^'^'-'^'■'-'^{M. x X) 

Sg is bounded in ^P^™ +i'^2'^"(M x X) and, as 5 ^ 0, converges to in 

^2m'+l+e-2l-n^^ X X) for all £ > 0. 



{SgU,U)g < 



{Egu, u) 



{Kgu, u) 



{Rgu,u)^ 



(u) is disjoint from the elliptic set of A'. This proves (8.6). 



Hence 
(8.7) 

so WF 

Finally consider (iii). We work by induction, assuming that the result is known 
for all nonnegative integers smaller than a given k. Note that Ay is an edge 
(pscudo) differential operator, so the conclusion of (iii) for positive k is only non- 
trivial if m > I + {n - l)/2. In particular p ^ WF™ ''(/iy) for m' sufficiently 
negative. We therefore prove the following analogue of (8.6) with u replaced by 
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If to' < ; + n/2 - 1, ue H-°^'\l X X) and p WF^' (m) then 



proceeding much as before. 

To prove (8.8), note that Lemma 1.6 gives a distributional equation for Uk = AyU 
of the form 



(8.9) 



Uuk + [QkD^ + ^Pfc) u = 0, where Pk G C°°([0, e); Diff2'=+i(y)), 

QfeeC~([0,e);Difr2'=-i(F)). 



Now applying the test operator Ag and pairing with gives an estimate similar 
to (8.7) with an extra term: 

(8.10) \\A'gUk\\l- {Ssuk,Uk)g < {Esuk,uk)g + {KsUk,Uk)^^ 

{R5Uk,Uk)g + (AsA'^u^AsBku)^ 

The first term on the left can be reorganized, modulo terms uniformly controlled 
by the inductive hypothesis, to bound a positive multiple of \\A'gu\\H2i<^{Y)' tiie tan- 
gential Sobolcv norm of A'gU. Wc now show that the last term on the right can be 
estimated by a small multiple of this norm, modulo the inductive bounds. Note 
that As = GsA'g + C's + Ds where Gs G ^'^^^^'^(R x X) is uniformly bounded as 
(5 J, 0, is lower order and Ds is supported in the region of known regularity. 
Modulo terms bounded by the inductive hypothesis or by the hypothesis of (8.8), 
we can thus write the last term as 



{A^A'su,BkG*sGsA'su) 



r2fc/ 



We may rewrite Bk as a sum of terms x~^CkSi with Ck G Diff {Y) and where 
Si are smooth b- vector fields. Hence modulo controllable terms the inner product 
above is estimated by a sum of terms of the form 

\\A^A'su\\\\x-^CkSiG*sGsA'su\\^. 

The latter norm is (again modulo known terms) controlled by e||^^M|j^2fc(Y), with 
the e coming from the small support in x. Thus the last term in (8.10) can be 
absorbed in the left, giving the inductive estimate and proving (8.8). 

To prove (iv), we appeal once again to Lemma 1.6 and proceed as with (iii). □ 

Theorem 8.3. Letu{t) be the solution of the Cauchy problem (4.1)-(4.3) and sup- 
pose that u e H~°°'\l X X) for some open I 3 i and some I G M then 

p e i?±,o(*). r'(p) n WF™(u) = 9 for m < l + {n - 1)/2 and e > O 

=^p^WF"'-\u) V(5>0. 

Proof. By (i) and (ii) of Theorem 8.1, there is no edge wavefront set of order 
{m,m — (n — l)/2 — S) along all bicharacteristics in ^S^^g-^i^ ^ ^) terminating 
at t = t, y = y{p)- Hence by (iii) of Theorem 8.1, closedness of edge wavefront set, 
and Hormander's theorem, p ^ WF'"~^^(u). □ 
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9. The rescaled FBI transformation 

To analyze the boundary behavior of solutions to the wave equation, we employ 
a rcscalcd version of the onc-dimcnsional FBI ("Fouricr-Bros-Iagolnitzcr") trans- 
formation in the time variable. The FBI transform is a variant of the Bargmann 
transform which was employed by Bros and lagolnitzer [10] to study microlocal reg- 
ularity in the analytic setting. Its properties were further elaborated by Sjostrand 
[21]. For our purposes, the FBI transform could be dispensed with in favor of a 
composition of localization in time and Fourier transform. Use of the FBI trans- 
form, however, seems more likely to admit generalization. The C°° approach to 
the FBI used here follows the spirit of [24], which deals with compact manifolds. 
A partial FBI transform is also used in a related analysis of the wave equation by 
Gerard and Lcbcau [6]. 

Consider the complex phase function ^(f, r, t') = i{t — t')^{T)/2 + {t — t')T. The 
associated FBI transform applied just to the time variable is 

Tu{t, r, x,y)= [ e"^^*'^'*'^a{t, t' , T)u{t' , x, y) dt' 
Jx 

where a is a polyhomogeneous symbol of order 1/4 with proper support in t, t' . It 
is an elementary consequence of the stationary phase lemma, demonstrated in [24] , 
that T*T is a pseudodifferential operator on K of order and hence is L^-boundcd. 
It is elliptic at ±oo) £ S*M. if a is elliptic at t = t' = t, t = ±00. Thus the 
cUipticity of a at {t, 00) implies that there exists G € *°(IR) with {t, 00) ^ WF'(G) 
such that 

(9.1) T*T = Id +G, so Tm = in T > 1 near {t = t} ^ u = -Gu near {t = t}. 
The basic intertwining property of T corresponds to the boundedness of 

(9.2) TDt - tT : ^^(K) — > (r)^i^(M^). 

We shall choose a to have support in r > 1, the region r < — 1 being handled by 

reflection in t. Fixing a product decomposition near the boundary and inserting a 
cutoff x(a;) with x(a;) = 1 for x sufficiently small and with support in the product 
neighborhood of a boundary component Y, we scale the a;-variable and define 

(9.3) Su{t, T, x,y)= [ x(-)e'*(*>->*')a(t, t', T)u{t', -,y) dt'. 

J T T 

If we let X denote the model cone [0, 00) x Y then 

(9.4) S :C-°°{R X X) — >C-°°(R x [L oc.) x X). 

In fact the assumed properness of the support in t, t' ensures that Su is C°° in r 
and t. However, it is the growth in r, uniformly in x, that will interest us here. In 
view of the assumption on the support of a is we will use a decomposition of the 
type described in (4.6) of m G C~°°(R x X) under convolution in t and examine u+. 

Taking into account the scaling, the boundedness of T implies that S is bounded 
on the scale-invariant spaces in the x variable 



(9.5) 



S : Lg(M X X) — > Ll{M. x [1, 00) x X) 
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where the measure on the left is dxdtdy/x and on the right dtdrdxdy/x and we 
assume that a is a function of t — t' only to get the global estimate. The scale- 
invariance of the measue shows that 

(9.6) S*S = x{x)G 

where G is the pseudodifferential operator on M in (9.1). 

The simplest intertwining properties of S follow directly from (9.3): 

SoQu = QoSu, Q€ Difr*(y) and 
S{xu) = —Su. 

T 

It follows from the second of these that if / € C°°{X) then composition with / as 
a multiplication operator gives 

(9.8) Sof:Ll{RxX) — ' f] (~) LI{R x [l,oo) x X). 

Note that x/t is bounded on the support of Su and r > 1 by assumption. Since 
we will only use the FBI transform to examine behavior in Taylor series at the 
boundary, we shall denote the space of operators with the property (9.8) as E°° 
and consider these as error terms. The presence of the localizing function x in the 

definition of S means for instance that 

(9.9) SoxD^ = xDioS + S'o /', /' G C°°(X), S' o f g E°° . 

Here, S' is an operator of the same form as S with a different cutoff x- 
Scaling (9.2) and using (9.7) we conclude that 

(9.10) S o {xDt) - xS : Ll{M. x X) — >t-^{x)LI{Rx [1,oo) x X). 

The scaling of x makes x a global variable on IR+ , and we arc primarily concerned 
with behavior as i; — * oo. The structure which arises at infinity here corresponds to 
the IR+-homogeneous metric in (1.13). This is a conic metric on X = [0, oo) x Y, but 
its uniform behavior near infinity is rather dififerent from its behavior near i = 0. 
Consider the inversion 

(9.11) X° = {Q,oo)xY 3 {i,y) i — > {l/i,y) G (0,oo) x F = X° 

which induces an isomorphism oi C~°°{X). Under this transformation, the metric 
becomes a scattering metric in the sense of [16] at the boundary w = 1/x = Q, 
which is to say it is a particular type of asymptotically locally Euclidean metric. 

In [16] the associated compactly-supported Sobolev spaces H™^{X) are defined 
for any manifold with boundary. It is natural then to introduce Sobolev spaces 
which arc of "b-type" near x = and of "scattering type" near .x = oo, based 
however on the R_|_-invariant measure dxdy/x (note that this convention differs 
from the weight used in [16]). Thus if we choose <p G C^([0,oo)), with (f){x) = 1 
near x = then we may define 

(9.12) H^.,,{X) = {ue C-°°(X); <^(^)« e H^^X), 

(1 - <P){l/w)uil/w,y) e H^,,iX)}, m e R. 

Since both the "b" and the "sc" Sobolev space reduce to the standard Sobolev 
spaces in the interior, this is independent of the choice of (p. We shall also employ 
weighted versions of these spaces which we will generally write in terms of the 
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weights X and {x). The former is a boundary defining function near but is also, 

near infinity, of the form of 1/w where ?i> = is a defining function for inverted 
infinity. On the other hand is just a defining function for inverted infinity. 

Thus we consider the weighted spaces 

(9.13) x\x)^H^_^^{X), l,k,m&^. 

Directly from the definition, these scales of Hilbertable spaces satisfy complex in- 
terpolation. 

The spaces which arise here correspond to functions of r and t with values in 
these weighted Sobolev spaces. We are primarily interested in global behavior in r 
but local behavior in t. We therefore define 

(9.14) LL(Mt X [1, oo]; H1^.^^{X)) = {u G C-°°(IR x M x l);u = in r < 1 and 

ct>{t)u eL\Rx R;H^_JX)) V G C-(M)}. 

More generally weights will be written out as in (9.13). Note that the closed bracket 
at r = oo in (9.14) is intended to indicate that these spaces are indeed global in r. 

Lemma 9.1. For any m, I E M. 

(9.15) S:x'h:^{RxX)^ f| rt'(i)-"LL(K* x [1, oo]; i?;?-"(X)). 

ae[0,m] 

Conversely provided x' € C°°(X), x'x = x' o,nd m > 0, 

(9.16) ifu&Ll{M.xX)then 

5« G fl r-'^'(i)-"iL(Kt X [1, oo]; i?,"-«(X)) near {t = t}^ 

ae[0,m] 

xV = ui +U2, ui G a;'iJ™(K x X), U2 G C~(M;L^(X)) near {t = t}. 

Proof. For positive integral m, u E .t'H™(IR x X) with support near the boundary 
if {xD^)''{xDt)PQu G x^Ll{RxX) for aU Q G Diff*(r) and all A;+p + ord(Q) < m. 
The continuity estimates in (9.15) then follow from (9.5), (9.7), (9.8), (9.9) and 
(9.10). Similarly, for negative integral to, u G H1P{M. x X) may be written as a 
finite superposition of edge operators of order —to applied to elements of L'^{R x X) 
and continuity follows similarly. Complex interpolation on both sides then gives 
the general case of (9.15). 

To see the partial converse, (9.16), first replace u by u+ as in (4.6). Since the 
amplitude of S is assumed to be supported in r > 1 this simply changes S by a 
similar operator with rapidly decreasing amplitude. Thus we may suppose that the 
condition holds for u+. The invertibility of S in (9.6) shows that S*S = x modulo a 
term arising from G; applying the arguments above to S* therefore gives (9.16). □ 

It follows from (9.15) that if G C~(M) is 1 near 

(9.17) {l-(f){x))SoB -.x^H^iRx X) — > 

T-\x)-^ H^^^{R X [l,oo) x X) V M, Bg *-~(K x X). 

Thus Su is rapidly decreasing in the sense of Schwartz near infinity. Thus the 

boundary part of the scattering wavefront set, in the sense of [16], of Su is deter- 
mined by the edge wavefront set, over the boundary, of u. In fact we need to use 
a notion of scattering wavefront set in a uniform sense in t and r. Such uniform 
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versions of the scattering wavefront set are also discussed in [16] . As there, we are 
only really interested in the scattering wavefront set near x = oo. 

On a compact manifold with boundary, X, the scattering wavefront set is a subset 
of the boundary of the radial compactification of the intrinsic scattering cotangent 
bundle. This bundle reduces to the usual cotangent bundle in the interior but has 
basis at any boundary point p the differentials dx/x^ and dyj/x, where where x is 
a local boundary defining function and y = at p, are boundary coordinates. The 
part of the scattering wavefront set over the boundary can be defined explicitly as 
follows. If (C 77) G ^^T*X is the coordinate representation of a general point in the 
fiber over P, Cfi + ^ ' ^ e """T^X then (C,r?) ^ WFsc(u) if for some V e C~(X) 
supported in the coordinate patch and with ^(p) 7^ 0, 

(9.18) j e-'i+'^^^ix, y)u{x, y) ^ dy is C°° near (C, r?). 

Note that changing to the singular variables l/x,y/x identifies the complement of 
the boundary in a neighborhood of p with an open set in M". Then the formal 
integral in (9.18) may be identified with the Fourier transform. Since u is an 
extendible (hence Schwartz) distribution, (9.18) is well defined; this condition is 
also independent of the choice of coordinates. 

The uniform version of the scattering wavefront we need is obtained by strength- 
ening (9.18). 

Definition 9.2. If u G T-'Lf„^(R x [1, 00]; (i)PH^^jX)) for some m, / and p then 

(9.19) (f, 00; y, C, n) i WF,e,a(t;; t-'L^J^ x [1, 00]; l,xfHl^{X))) 
\iv = v\+V2 where, for some m, v\ e t~^L\^J^ x [1, 00]; (xY Hl^^jJC)) and 
(9.20) 

/I r/r 
^^^rfy g X [1, 00]; C-(M")) near {t,C,fj), 

x X 

where ip{x,y) G C^([0, 00) x Y) is equal to 1 near (0,^). 

Note that the regularity, m, here is irrelevant since it is the boundary part of the 

scattering wavefront set wc arc examining. 

We will show that the wavefront relation of S is related to 

(9.21) 7^ = {(C',2/,M,^;e,2/^r/,^^A);A>0, t = t', y = y\ (C, m) = A-i(^, -r/)}. 
Proposition 9.3. If u & x^H-^{U. x X) has WF~''(u) c {A >} then 

(9.22) W¥,,^a{Su; t-'lI,{R x [1, 00]; (x)'-"" Hl,,{X))) 

c n o (wftHu) n ^5^xax(K X X)) , 

in the sense of Definition 9.2. 

Proof. If u e x''H™{R X X) then, by Lemma 9.1 the scattering wavefront set of Su 
relative to T~^Lfg^{R x [1, 00]; (i)'~™if^.gp(X)) is empty, since Su lies in the space 
T-^Ll^{R x [1,00]; (.i)'""i7{;_,^(X)). Thus it suffices to prove (9.22) for to = 00. 

It is therefore enough to consider u G x'ffg" (M x X), for some to, having support 
in the product neighborhood of the boundary and having scattering wavefront set 
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contained in {A > 0}. Then, using a partition of unity in the edge calculus, we may 
decompose 

u = u' + BjU 

3 

where u' £ x'^H^^ and the Bj have small wavefront sets, in A > 0. Again by 
Lemma 9.1 the term u' produces a term which is Schwartz near x = c» as a function 
with values in the weighted space t^' Lf^^{M. x [l, oo)). This corresponds to the 
absence of any scattering wavefront set. So we may replace u by Bu, with an edge 
pseudodifferential operator, B e vE'°(K x X), with essential support concentrated 
near some boundary point (t, 0, y. A, ^, ^). We may also suppose that B has its 
support near the boundary component Y. Since the edge wavefront set is a conic 
notion this means that we may take A = 1. In local coordinates 

(9.23) Bu{t, X, y) = (27r)-i j e^((«-i)«+(y-y')'/+TA) 

b{t, X, y, s, y', T, ^, r/, X)u{xs, y', t — xT)x^ ds dy' dT dr] dX, 

where, at the expense of another error of order — oo in the edge calculus, we may 

assiime that 

the amplitude b has small conic support (as A — > oo) around 

(9-24) - _ , _ 

t = t, x = 0,y = y, s = 1, y =y, T = 0, ^ = and r] = rjX. 

Since we are interested in the scattering wavefront set of the image we can localize 

near x = oo; the assumption on b means the image is already effectively localized 
in y near y and in t near t. Thus, in local coordinates based at y = 0, we take the 
(normalized) Fourier transform of v = {1 — (f){x))SBu, 

(9.25) J'v{$,',ii;t)= j e'^^'+^y^'v{x,y;T)x''-^dxdy. 

We then localize near some point using a cutoff <f>{x,y, , ii) to examine the reg- 
ularity; such regularity allows us to conclude absence of scattering wavefront set 
at points in the support of (p, where ^' and /i are coordinates in the fibers of the 
scattering cotangent bundle determined by the canonical one-form 

^'dx + fi ■ d{xy). 

Finally, then, we arrive at the composite operator 

(9.26) Gu = T{1 - (j){x))SBu, 

/TS tT 
e'^cu{—,y', t' )T-2j"+ids dy' dT dt' dx dy d^ drj dX, where 
T T 

c = c{x,y,^' , l^,t,t' ,T, -, s,T,^,r], X) 

T 

~ X 

= Hx, y, tJ)a{t, t', T)b{-, y, t', s, T, ^, r], A) and 

T 

i; = £e + xyfi+{t- t')T + ^{T){t - t'f + {s- l)e + {y-y')-ii + TX. 

First suppose that e{x, A) is a conic cutoff keeping \x—X\ < 6{x, A) asymptotically 
for some small 5 > and equal to one on a smaller region of this type. Inserting 
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1 — e{x,X) into the integral we may integrate by parts using the vector field 

(9.27) Vi^dT + -df satisfying Vii^ = X - £ + —^{t' - t). 

T T 

Since t — t' is small, for large r this is elliptic in A and x on the support of the 
amplitude. This gives an operator of the form (9.26) but with amplitude arbitrarily 
decreasing in x, A. Rapid decrease in A implies rapid decrease in rj and ^ and hence 
again effectively replaces B by an operator of arbitrarily low order. Hence this term 
also has no uniform scattering wavefront set. 

Thus, by inserting the cutoff e wc may assume that a; ^ A on the support of the 
amplitude. Now, consider the vector fields 

(9.28) V2=ds- -ch - —Tdv and W = dy 

' 3 TS 

which annihilate u{^,y',t' — ^) and satisfy 

(9.29) = <e - — - + i^{t' - t)) and Wi; = rj + xfx. 

S S S \ T J 

These functions, with t — t, and Vitp, define TZ. 

Assuming that the support of (f) in (9.26) is disjoint from the image of the essential 
support of B under TZ, integration by parts now allows the amplitude, c, in (9.26) 
to be replaced by a symbol of arbitrarily low order in ^, rj and A. Simpler versions of 
the arguments used to prove (9.15) and the boundedness of edge pseudodifferential 
operators show that this leads to a map 

(9.30) x^H^ci^ xX)^ C°°(IR" X Mj; t-'l'^{[1, do)). 

This proves (9.22). □ 

We shall apply the FBI transform to solutions of Du = 0. Since S" is a Fourier 
integral operators, with complex phase function, the transform satisfies a model 
equation involving the Laplacian, Aq, on the model cone X = R_|_ x Y with respect 
to the the product metric dx^ + x'^hg^y, dy). 

Lemma 9.4. There are operators Si and S2 of the same form as S with amplitudes 
of order j and L e Diff^(X) such that 

(9.31) {Ao - 1)S = -T-'^Sn + T-iSi + T-^x-'^SiL. 
Proof. Computing directly 

ntjju = ipnu - 2{D^tjj){D^u) - {Dlil))u. 

Let S be defined as S was but with cutoff G C°°(X), still supported in the product 
neighborhood of the boundary, with — ip. Then 

(9.32) 5(n^w) = S{au) + S{Liu) with Li e DiSl{X). 
Now the commutation relation for Dt shows that 

SD^iju = T'^Stlju + T^S'u 
where S' is of the same form as S but with a different amplitude of order j . Since 
SAo = AqS, AtP - Aoip e a;-i Diff^(X) and S{tpu) = Su we conclude that 

(1 - Ao)S{u) = T-^SUiju - T-iS'u + S{A - Ao)ipu, 
yielding (9.31). □ 
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10. Reduced normal operator 

The reduced normal operator of the conic wave operator, as an element of the 
weighted edge calculus, is 1 — Aq where Aq is the Laplacian for the tangent product- 
conic metric on X = [0, oo)j x y. Regularity results obtained in subsequent sections 
depend on global invertibility properties for this operator. It is important to un- 
derstand that its behavior at infinity in the product cone is quite different from 
its behavior near i = 0. In fact, near oo, the metric (1.13) is a special case (again 
the model case) of a scattering metric considered in [16]. The local and microlocal 
estimates in [16] arc combined below with the analysis of the domain for conic met- 
rics above to get the desired invertibility estimates, which are of a standard type 
in scattering theory. 

The analysis of the Fricdrichs extension of Aq on X proceeds very much as 
above, using also the fact that the metric is complete at infinity. Thus the scattering 
Sobolev spaces are, near infinity, the natural metric Sobolev spaces and the domains 
of powers of Aq are closely related to the "b-sc" Sobolev spaces defined in (9.12). 

Proposition 10.1. For the Friedrichs extension of Aq, the Laplacian for the met- 
ric (1.13), 

Dom(ZiJ) = x-^+P{x)-^HlJX), <P < ^, 
(10.1) ^ ^ ^ 

Dom(Zio*+') = x^{x)-^-^Hl^^\x) + Ccf>,0<S<So 

for Sq > sufficiently small and with (j) G C^(X) identically equal to 1 near x = 0. 

The spectrum of Aq is the whole of [0, oo). Just as in (3.23) in the compact-conic 
case, Aq defines a continuous map, 

(10.2) 

Aq : x^{x)-^-^Hl^f{X) + C(/) ^ x-^+^{xf-^-^Hl~^+\x), S>0 small , 

however this is never an isomorphism. On the other hand, the spectral family 
Z\o — A, for A e C \ [0, oo) does define a continuous map as in (10.2) which is an 
isomorphism. 

One of the standard results of scattering theory, proved in [16] in the context of 
scattering metrics, is the limiting absorption principle. This asserts the existence of 
the limit of the operator {Aq — X±ij)~-^ for A G (0, oo) as 7 | 0. Of course the limit 
cannot exist as a bounded operator inverting the resolvent, i.e. cannot be defined 
on the range space of (10.2). However it does exist on a somewhat smaller space. 

Proposition 10.2. For 7 > the resolvent 

(10.3) (Zio - 1 ± ^7)-' : £-'+'{if-'-^Hi;^'+\x) ^ 

x'{x}-'-^Hl^{X)+CcP 

restricts to an operator 

(10.4) {Aq-1± i^r' : x-'+'{xf-'-^-'Hl-'+\x) — 

x'{x)-'-^-^Hl^^\x)+C<f> 
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for any s and as an operator 

(10.5) (Z\o - 1 ± 11)-' : i-^+\if-'-'i--^-'Hi:^^\x) 

~5^~p5-f +i+-ijt+^(X) + C0, for any e > 0, 

the strong limit exists, os 7 J, 0. 

Proof. The only difficulty with the convergence of the resolvent is related to the 
large, i.e. scattering, end of the cone. Modulo compact errors the analysis in [16] 
therefore applies and gives the same continuity properties for the limit of the resol- 
vent as in the scattering case itself. □ 

Thus we have two limiting operators (Z\o — 1 ±«0)^^ with the mapping property 
(10.5). For these operators it is not possible to take e < in the domain space, nor 
in the range space without further restriction. However there are "larger" spaces 
on which the limiting resolvent is defined; these are fixed in terms of conditions on 

the scattering wave front set. 

Proposition 10.3. Let ^' he the dual variable to dx in the scattering cotangent 
bundle. Let U± be an open neighborhood of the part of the radial set for Aq — 1 on 
which > 0. Then for any m e R and k' > k > 

(10.6) 

{Ao-l± iO)-' : {/ e i-2+^(i)'-*-*+'=if|;/+'(X); WF,e(n) n C/± = 0} ^ 
{u e x^{x)-^-^+'''hI+\x) + C<^; WF,e(w) n C/± = 0} 
and this is a two sided inverse 0/ Z\o — 1 on the union over k and k' of these spaces. 

Proof. This just follows by combining the results of [16] with the analysis above of 
the domain of Aq. All the estimates of [16] concern the behavior as x — > 00 and 
the appearance of the conic boundary at a; = makes essentially no difference to 
the argument. This result could also be proved using the analysis of the forward 
fundamental solution of the M+-invariant conic wave operator from [2]. □ 

11. DiFFRACTIVE THEOREM 

Proposition 11.1. Suppose u G C(]R, P^-^), for some 6 > 0, is an admissible 
solution to the conic wave equation and 

(11.1) WF(u) n i?|(0) = 

for some e > 0. Then for any 4> G C^(M x X) with support in a sufficiently small 
neighborhood of {0} x dX 

(11.2) <pu e x-^'H^{m. X X), y s' >o, 

so u is conormal near the boundary for \t\ small. 

Remark 11.2. Note that, since 21 is a solution of the wave equation, (11.2) is equiv- 
alent to the condition that for all k, (p{x)D^u is an function of t near t — with 
values in f]^Vs, provided (j) G C°"{X) is chosen to have support sufficiently near 
the boundary, i.e. the conclusion of Theorem 4.4 holds with s = 00. 
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Proof. First we replace u by u+ in terms of the decomposition (4.6), this has the 
same regularity properties but its Fourier transform in t has support in r > 1. The 
treatement of U- is similar. Now we proceed in three steps. 

First, the hypothesis (11.1) allows Theorem 8.1 to be applied directly. Since 
u e C(M,£n) c a;-*Lg_iojM x X), the first part of Theorem 8.1 shows that 

(11.3) IC(0) n WFr'-'^(w) = V m, V ^ > 0. 

Our initial goal is to apply the identity in Lemma 9.4 to show that 

(11.4) (t>D'lu&x-^H-°°{^x X) 
for all k and a cutoff as in the statement above. 

X ——6 

The initial hypothesis on u implies in particular that u G x~°Hi (R x X) near 
{0} X dX. Applying Lemma 9.1, with the cutoff in the definition of S supported 

sufficiently near a; = we therefore conclude that 

(11.5) 5ue fl i-^(i;)-VLf„,([l,oo)^xRt;ir|;/""(X)) nearf=0. 

c«6[0,n/2-<5] 

Similar conclusions apply to the terms on the right in (9.31), with S\u having the 
regularity property (11.5) and 

(11.6) 

S^{Lu) g fl ([l,oo). X M,;ff|;/-'-"(X)) 

a6[0,ra/2-<5-2] 

near f = 0. 

It is also the case that near i = 0, Su is in the domain of Z\q if (5 > is small 
enough. The identity (9.31) then holds, near x = 0, in the sense of the domain of 
Aq , which is just a weighted b-Sobolcv space. 

Of fundamental importance is the estimate on the scattering wavefront set which 
follows from (11.3), together with the initial assumption (11.1), by use of Proposi- 
tion 9.3. Namely, the incoming wavefront set for all terms, Su, Siu and 5*2 (iu), 
computed with respect to the spaces in (11.5) and (11.6) is absent — these func- 
tions are microlocally rapidly decreasing as a; — > oo in this uniform sense. Now, 
Lemma 10.3 shows that solving (9.31) gives an "improvement" in the estimate 
(11.5) by a factor of t~2 at the expense of more growth as i oo. Since we may 
invert S, up to errors which involve a smoothing operator int, we may iterate this 
argument, reapplying Theorem 8.1 to obtain absence of scattering wavefront set 
with the new t weight and using the improved estimate on Su to estimate Siu and 
S2{Lu). Such iteration yields an estimate 

(11.7) 

Su e f]x-^{x)^T-^Hl~\R X [1, oo) X X) near F = 0, VM with N = N{M). 

OL 

Using the reverse regularity estimate in Lemma 9.1 we conclude that (11.4) does 
indeed hold. 

Finally we apply Theorem 8.1 again to deduce the claimed regularity (11-1). 
Indeed, we may apply the first three parts of the theorem, in succession, to D\u 
for any integer k. As before we find that over t = and the boundary 

(11.8) WFf '-*(u) c OG(0), V M. 



44 



RICHARD MELROSE AND JARED WUNSCH 



Now we may use the outgoing propagation result to conclude that 

OG(0) n wf^-^'-\d^u) = 

for any 5 > 0. Since this is fixed finite regularity, certainly implying that D^u G 
Ll{R X X) locally for all t derivatives, (11.2) follows. □ 

Proof of Theorem 4-4- Wc use the decomposition (3.18) to assume without loss of 
generality that there are no singularities in {t < 0}. By time-translation invariance, 
we may assume that i = 0. Choose a cutoff tp{t) G C°°{R) which is 1 in t > —^5 
and in t < —6 for 6 > which will be chosen small. Then v = il){t)u satisfies 

(11.9) nv = f = -2i'4)'{t)DtU - tp"{t)u. 

Now, if (5 > is small enough, the hypothesis on the wavefront set of u means that 
WF'^~^(/) is disjoint from all incoming rays arriving at the boundary at time 0. 
Thus, it may be divided into two pieces / = /i + /2, with both /j supported in 
— j5>t>—25,fi supported away from the boundary and satisfying 

(11.10) /i e Cf (R; Vr-i-p), yp&n 

and with WF(/2) disjoint from the incoming cone at t = 0. It follows that the 
forward solution to Ovi = /i is in H^(M. x X) away from the boundary and is in 
C{M.; Vj.) near it. Thus it suffices to consider — /2 which is a solution near t = 
with no incoming singularities at all at time t = 0; by the uniqueness of solutions 
it is equal to u — Vi near t = 0. Now V2 may be extended to a global solution u' 
which is equal to W2, and hence to u — vi near t = and this too has no incoming 
singularities at time 0. After a finite amount of smoothing in t using Lemma 4.3, 
Proposition 11.1 applies to u', showing that it has no outgoing singularities at time 
(and indeed lies in C{R;T>oo)) locally. □ 

12. Propagation of tangential regularity 

As a prelude to the division theorem, we will prove that regularity of solutions 

in the tangential, i.e. boimdary, directions is conserved imdcr time-evolution. This 
result represents a global version, in energy spaces, of the microlocal estimates in 
Theorem 8.1, part iii. 

Recall that for a boundary component Y oi X, Ay denotes the Laplace-Beltrami 
operator induced on Y by the metric /iq. Set 

Ys = (I + AyY^^, seM. 

These tangential pseudodifferential operators act naturally on the boundary, but 
may be viewed as acting on the fibers of the product decomposition in Theorem 1.2. 
Thus they act on functions or distributions on (a, b)t x dX n [0, e)x- 

Recall that denotes the energy space Vg © I?s-i where V = Dom(Z\^/^). We 
write £ = £i as well as V = Vi for convenience. £i is a Hilbert space with the 
norm 

\\{u,v)\\l = \\u\\l^ + \\du\\l^ + \\v\\l^, {u,v)G£. 
Note that if we let 

denote the infinitesimal generator of U{t), M is not selfadjoint with respect to the 
norm on £ as defined here (owing to the ||u||^2 term). 
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For use in the sequel, we note the following energy estimate: For an admissible 
solution (u, Dtu) with {u, Dfu) G C°°(M; £00), and an operator Q such that 

(12.1) Q:C°°(M;Poo)^C°°(M;P„/2) 

set 

Qu{t) 



Dt{Qu{t)) 



We can then compute 



iiii*(oiii^ae((:^g:)),. 



£ 



(12-2) =Re(^M*,*), + Re((^py 

= Re {i[n, Q]u, A(Qw))z,2 + Re {iDt{Qu), Qu)^^^. 

Let denote a product neighborhood of the boundary and let £c{0) C £ and 
PcC^) C f be the subspaces with compact supports in O. 

Lemma 12.1. // B{x) is a tangential pseudodifferential operator of order 1 de- 
pending smoothly on x and B{0)C = then x~^B{x) : VdO) — > L^{X). 

Proof. Wc know that operators of the form xr^V with V a tangential vector field 
are bounded in this way, as is B{x) itself. The generalized inverse of the tangential 
Laplacian satisfies EAy ~ Id— Hq where Ho is orthogonal projection onto the 
constants. We have B{x)IiQ = xB{x), hence 

(12.3) B{x) = B{x)EAy + B{x)nQ = ^Bj{x)Vj + xB{x) 

3 

where Vj are vector fields tangent to dX, the Bj are of order (and hence bounded 
on L^) and B{x) is again of order 1 and smooth in x. Thus 

(12.4) x-'^B{x) ^^Bj{x)x-^Vj + i3{x) 

3 

is bounded as claimed by (3.9). □ 

Let C C O be a smaller product neighborhood of dX. 
Proposition 12.2. For all \t\ sufficiently small and s e M, 

YsU{t)Y-s ■ £c{0') £ 

is bounded. 

Proof. As {C°°{X)+C)®C°°{X) c £00 is dense in £, it suffices to prove the relevant 
estimate for Cauchy data in £ao- 

Let $(t) = {u{t) , Dtu{t)) = U{t){uo,ui) be the solution to the Cauchy problem 
(4.1)-(4.2) with supp WO) suppui C O'. Then there is an open interval I containing 
such that for t G /, wc have suppu(f), suppDtu{t) C O. The hypothesis (12.1) is 
satisfied for Q = Yg, hence 

(12.5) ^ ^ WYMt) \\l = Re {Y,ut, Y,u) + Re {[A, Y,]u, Y,ut) 
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for all t & I. By Lemma 1.7 and Lemma 12.1, ||[/4, VsJuH is bounded by a multiple 
of ||(w, 0)||£., locally near Y. Hence 

~\\Ysm\\l<c\\YMmi 

from which the boundedness follows. □ 



13. Global propagation of conormality 

The tangential regularity discussed in the previous section is the main step to 
showing that incoming conormal waves with respect to the surface R±j propagate 
through the boundary to be conormal on the outgoing radial surface R±,o- To 
prove this we need further to show that regularity is preserved under the repeated 
action of the radial vector field 

R = xD^ + {t- t)Di. 

where t = t is the time at which the cone hits the boundary. Using time-translation 
invariance we may always assume t = 0. 

To begin, note that we can combine identities from Lemmas 1.6 and 1.7 to obtain, 
for any fixed s,k, 

(13.1) 

fe-l k-l 

[□, YsR''] = J2 CjYs+iR^D + + x-^Pj)Y,+iW + {RD^ + x-^S)Y,R'' 

j=0 j=0 

where Cj G C, Oj G C°°([0,e) x Y) and R e C°°([0,e) x Y^^^-^Y)), and where 
S, Pj G C°°([0, e) X Y;'^^{Y)) annihilate constants at a; = 0. Since the opera- 
tor YsRk satisfies the hypothesis (12.1), (13.1) makes sense when applied to an 
element of C°°(M; 2? oo), with the equality holding in C°°(M; I?„/2-2)- Thus using 
Lemma 12.1, the energy identity gives 

(13.2) ~\\{YiR^u,DtYR^u)\\l 

<C\\{YiR^u,DtYiWu)\\^ \\{Yi,R^'u,DtYi,R^'u) 

i'+j'<i+j 

Proposition 13.1. IfU{t) is the solution operator to the Cauchy problem (4.1)- 
(4.3) then for a product boundary neighborhood O, small time T and each fc € N, 

(13.3) 

^ \\YiR^Uit){uo,ui)\\£<C \\YiW{uo,ui)\\£, \t\<T, iuo,ui) e £^{0). 

Proof. We apply (13.2) inductively. Note that 

Y \\YiWU{t){u{t),Dtum\s-^ Y \\{yiR'u{t),DtYiWu{t))\\s 

i+j<k i+j<k 

since 

i-i i-i 
[A, YiR^] = Y ciY^R'Dt = YiY c'lDtR'. □ 
;=o ;=o 
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Proof of Theorem 4-8. By results going back at least to Hadamard (see also [8, 4]), 

under the hypotheses of the theorem, there exists e > such that for < t < e, u{t) 
and Dtu{t) are conormal distributions in R x X at {x = a; — U {x = a; + with 
respect to H^{R x X) and H^~^{R x X) respectively. The distribution conormal 
to {x = X + 1} remains conormal for small positive time, as it does not reach dX; 
we thus assume without loss of generality that this component vanishes. Hence for 
all s, k 

Y,R\u,Dtu)eC{Io;£s) 
for some interval Iq containing but not necessarily x. 

Note that with from Definition 4.2, [8s,-R] is a properly supported pseudo- 
differential operator of order s on R, hence 

(13.4) [es,R]=Qes+E 

where Q E ^'''(R) and E G ^'~°°(M). Since convolution with a properly supported 
time-translation invariant pseudodifferential operator of order on Rt maps finite 
energy solutions to finite energy solutions, iterative application (13.4) shows that 
for all i and j 

where 7 C /q is a time interval containing 0. 

Now by Proposition 13.1, YiR^n^Qs-i{u,Diu) G £ for aU fc G N and for 
all t < T. The symbols of Yi, R, and □ are defining functions for the conormal 
bundle to the hypersurface x = \t — x\. Hence {<ds-iu{t),Qs-iDtu{t)) is conormal 
to X = \t — x\ in (H^, L^). By Lemma 4.3. {u{t), Diu(t)) is conormal to x = |< — x| 
in (i?*, ff*~^). The theorem then follows by restriction to fixed t. □ 

14. Proof of the division theorem 

In this section, we prove Theorem 4.7. We begin with a preparation argument, 
allowing us to replace our hypothesis of regularity at R^ j with global regularity 
of the same order. As in the proof of Theorem 4.4, we may assume without loss of 
generality that u has edge wavefront set only in {t > 0}. 

Let z : X X° ^ R x X° be the inclusion map, and i* the induced contravariant 
map on cotangent bundles. Under the assumption of the nonfocusing condition 
(from Definition 4.5), there exists a microlocal neighborhood U of i*R'j{t) and 
k gN such that WF'^'^''{Y-i~u\t=o) fl C/ = 0. We construct a microlocalizer in such 
a neighborhood which preserves tangential regularity. 

Lemma 14.1. Let U C T*X° he an open conic neighborhood of i*Ri{t). There 
exist a smaller neighborhood V and an operator H G ^''^(X°), with Schwartz kernel 
compactly supported in an arbitrarily small neighborhood of the diagonal, such that 
WF' HcU, WF'(7 - 77) c V^, and [Yk, 77] = for all k e Z. 

Proof Let Ve e C^i^.) equal 1 for x < e/2 and for x > e. Let x e C°°(R) vanish 
for X < and equal 1 for x > 1. If e is sufficiently small then 

supp^,(|x - t\)M\v\Vm^ + C U; 

here rj is dual to dy and ^ to dx in T*X°. Let 

H = - t\)M^Y/{{Dxf + Ay))M\^ - 

where the function of the operator Ay/{{Dx) + Ay) may be constructed on the 
manifold x X and its Schwartz kernel cut off in x and glued into a product 
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neighborhood of {x = t}. Then H commutes with Ay, and hence with Yk for all 
k. □ 

Since the kernel of H vanishes near the boundary it preserves the domains of all 
powers of A. Let {v{t),Dv{t)) be the solution to the Cauchy problem (4.1)-(4.3) 
with initial data {Hu{0), HDtu{0)). By construction, WF(u - v) n Ri{t) = 0, so 
by the diffractive theorem, li — w G x~^Hj^{I x U) for some open interval I 3 t, 
U D dX. Hence it suffices to prove the desired results for v instead of u. 

Since 

Y_kiviO),Dtv{0)) = HY^k{u{0),DMO)), 
if the microlocalizer H is chosen concentrated sufficiently near the incoming set, 
the nonfocusing condition yields 

Y_k{v{0),Dtv{0))e£r+i. 
Proposition 12.2 now yields 

Y-kQs-iU{t){v{0),Dtv{0)) € £, 

hence 

Y^kU{t){v{0),Dtv{0))eCiR;£s) 
and the first part of the division theorem (that under the assumption of the nonfo- 
cusing condition) now follows by Proposition 5.5. 

We now prove the second part of the theorem (under the assumption of the 
conormal nonfocusing condition) . For simplicity, we now translate the time variable 
so that t = 0. 

We have 

RY_k{v,Dtv)\t^n 

= (xD^HY^kuiO) + tY^kDtuiO), xD,HY^kDiu{0) + tAHY^kuiO)). 

Since [H. A] and [H, xD^] are pseudo differential operators of order 1 and 0, com- 
pactly supported in X° , this is in £r+i provided the conormal nonfocusing condition 
holds. 

Now we prove the second part of the division theorem for the special case r+l = 1. 
By Proposition 13.1 

DRY_kV = f GCiR;V_i), 

hence 

{RY_kv{t), DtRY_kv{t)) = U{t)[{RY_kV, DtRY_kv)\t=o] + f U{t - s)(0, f)ds. 

Jo 

The first term in the right is in C{M.; £i) by the conormality assumption. The second 
is mC'^{R;£o), hence inC(E;£:i). Thus, RY_kv{t) e C(R;I?i), i.e. 

DtY^kV e t-^{CiR;-Di) + xD^C{R;Vi)]. 

By Lemma 3.2, we obtain for all e > 0, 

DtY.kV € t-iC(R;x-"/2+i-^L2(X)) C t-^^^-'x-''^^+^-'' LI{I x X) 

for all e > 0. We also know a priori that 

DtY.kV e C(R;Vo) C t^/^x-^'/^Ll.^^iR x X). 

Hence by interpolation, 

DtY_kV e a;-"/2+^/^-'Lg(/ x X) for all e > 0, 
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proving the theorem in the special case r + 1 = 1. 

If r + Z < 1, we apply Lemma 4.3 and (13.4) to conclude 

(14.1) Qr+i-iDtY_kV e a;-"/2+V2-^i2(7 ^ x). 

Thus 

DtY-kV e x-''/^+'/^-'H^+'-\l;Ll{X)), 
and since r + I — 1 < 0, Proposition 5.4 yields 

Dtv ex-''/'^+^/'^-'Hl+^-''-\l X X). □ 

15. Consequences of the division theorem 

In this section, we record two consequences of the division theorem. First, we 
prove a slightly weakened version of Theorem 4.6, the geometric propagation theo- 
rem. Then we prove Theorem 4.9, establishing conormality of the diffracted front. 

The weakened version of Theorem 4.6 we now establish is as follows 

Partial Theorem 4.6. If u is an admissible solution to the conic wave equoMon 
satisfying the nonfocusing condition in Definition 4-5 at {t} x Y for some r G M 
and I e (0, n/2) and if p G ^(f, Y), then 

(15.1) r (p) n WF^^'in) = 9 forr' €{r,r + l- 1/2) =^ p ^ WF'''-\u) V5 > 0. 

Ifu satisfies the conormal nonfocusing condition at {t} x Y, with k = l, then (15.1) 
holds for all r' G {r,r + 1). 

In other words, the theorem holds almost as originally stated, subject to limits 
on / and to the stronger conormal version of the nonfocusing condition, or yields 
one half derivative less, subject to the original hypothesis (again for a limited range 
of I). 

Proof. Let 

{v{t),Dtv{t)) = Qr+n/2-l{u{t),Dtu{t)). 

By (13.4), the new solution v satisfies the nonfocusing condition, with the same /. 
It also satisfies the conormal nonfocusing condition if u does, with r -\-l < 1. The 
division theorem now implies that there exists an interval I B i, and neighborhood 

D dX in X , such that 

Dtv e a;-"+'''-''+i/27j-oo(j ^ (j^^ 

where we allow r' e (r, r + / — 1/2) resp. {r,r + I) in the case where u satisfies the 
nonfocusing condition resp. conormal nonfocusing condition. On the other hand, 

the incoming regularity assumption shows that T^{p) n WF""/^^*" ~^{Dtv) = 0. 

Theorem 8.3 now yields p i WF""/2+'''"''"''(Aw), hence p ^ WF'''-''(w). □ 

We now prove Theorem 4.9 by including the effect of iterated regularity with 
respect to the vector field R defined in (1.17). 

Proof of Theorem 4-9. We will show that away from the geometrically propagated 
rays, the solution maintains its regularity under application YiR^ for all i e 2N, 

1 E N. The proof thus amounts to a microlocalizcd version of the argument used 
previously to prove conservation of radial conormality. 

Using Lemma 4.3, we may assume that r = 1. 
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By Proposition 13.1, for all i e 2N, j e N, there exists A'' = N{i,j) < such 
that 

YNYiR^ {u, Dtu) e C{R:£) 
for all t small. Thus, YiRWtU G x-^J?-°°(M x X), locally near t = t. 

We may apply the first two parts of Theorem 8.1 to the solution Dtu and con- 
clude that, for all m and e > 0, WF^'~'^{Dtu) is disjoint from all incoming bichar- 
acteristic segments, into and within the boundary, with endpoints p. Using the 
interpolation of wavcfront sets in Proposition 5.3 we conclude that the same holds 
for WF^ -'{YiRWtu) for all i, j and m. The third and fourth parts of Theo- 
rem 8.1. Then shows that p ^ WF e'"^''{YiR^u) provided m < -e. This establishes 
conormality. □ 

16. CONORMAL CaUCHY DATA AND THE PROOF OF THEOREM 4.6 

The most significant class of examples to which the conormal nonfocusing con- 
dition applies is given by solutions u with Cauchy data conormal to a hypersurface 
W which is at most simply tangent to the radial surfaces. 

Lemma 16.1. Let W he a compact smooth hypersurface in X° such that, in a 
collar neighborhood of the boundary, dx restricted to W vanishes at only finitely 
many points pj = (xj,yj), j = 1,...,M, at each of which W is simply tangent 
to X = Xj and suppose uq G H^{X) is conormal with respect to W then for any 
N gN, there exists k gN such that 

(16.1) Y_kUo e HUX°), . . . VpY_kUo G H^X"), p<N, s<r, 

for all Vi G C^{X°;TX) that are tangent to {x = xj} for all j. 

Proof Without loss of generality, we may localize and take M = 1. Since Y^k G 
\^-k(^Y), for any K we may choose k so that n{Y-k){y,y'), the Schwartz kernel of 
YLfe, is in C^(F x Y); we choose K>{n- l)/2 + N. 
The distribution u is of the form 



u{x,y) = j c{x,y,Oe"t'^-^y^U^, 



where c G 6'-''-V2+«(M x X°) (for all e > 0) is a symbol with one phase variable 
and ^ is a defining function for W. Hence 

Y_kU = j c{x,y',0'^iy-k){y,y')e"^^^''''^^ d^ dy'. 

U ^ ^ 01 dy(j) ^ on suppc, integration by parts using (p^^D^ or {£,<j>y>)^^ Dy> 
shows that Y-kU G H^^^{X°). Hence may assume that, on the support of c, there 
is one point, which in local coordinates we may take to be x = x, y = 0, at which 
(j) = Q = dy(j). The hypothesis of simple tangency and the Morse Lemma allow us 
to arrange that, locally, (j){x, y) = {x — x) — Y^- ajyj with aj = ±. The method of 
stationary phase now yields 

Y_kU = J (0-("-'^/'c(x,y,Oe'(--^)«d^ 

where c satisfies any finite number of the symbol estimates of any order greater than 
— r — 1/2 provided that K is taken large enough. Thus Y^ku e i7;+("-i)/2-^(^o) 
for any e > 0. Furthermore, as long as a+ |/?| < N, integration by parts shows that 
{{x — x-i)Dx)°'D^Y-kU remains in this space. □ 
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Such conormality for the initial data leads to solutions satisfying the conormal 
nonfocusing condition. 

Proposition 16.2. Suppose that {uq.ui) G £r{X) vanishes near dX and uq and 
Ui are conormal with respect to a hypersurface W C X° as in Lemma 16.1 then, for 
small f e M, the solution u to the Cauchy problem problem (4.1) -(4.2) satisfies the 
conormal nonfocusing condition with background regularity r and relative regularity 
I < (n - l)/2 to all orders k eN. 

Proof. By localization of the initial data we may again assume that there is just one 
point of tangency of W to one radial surface, x = c. Since the conormal nonfocusing 
condition is trivial microlocally away from the radial directions wc may assume 
that Mo and Ui are supported close to the point of tangency. Initial data for the 
wave equation which is conormal with respect to a hypersurface, W, gives rise, by 
Huygcn's principle to a solution which is conormal, for t > small, to the union of 
the two characteristic hypersurfaces W± through W; this follows from the original 
construction of Lax or from [8] . Here W± are each tangent to the radial surface 
{x ± t = c}. Thus W- is outgoing: the bicharactcristics forming it do not hit 
the boundary small times. Now, we may apply Lemma 16.1 above, regarding t 
as a parameter, to conclude that the solution satisfies the conormal nonfocusing 
condition with background regularity r for any I < (n — l)/2 to all orders. □ 

Proof of Theorem 4-10. We can write the fundamental solution in the form 



where m G X° lies sufficiently close to dX. Then for all e > 0, the initial data 



hence the solution lies is S-n/2+i-e for all e > 0. 

For small t > {Erh{t), DtEfn{t)) are conormal with respect to the hypersurface 
W which is the geodesic sphere of radius t around m. This is tangent only to 
the two radial surfaces x = d{m, dX) ± t and the tangency is certainly simple. 
Hence Proposition 16.2 applies and shows that u satisfies the conormal nonfocusing 
condition. Now Theorem 4.4 and Partial Theorem 4.6 shows that for, 2d(fri, dX) > 
t > d(rh, dX), the inclusion (L3) holds. On the diffracted front, but away from the 
direct front, application of Theorem 4.9 shows that E.^ is conormal and of Sobolev 
regularity — n/2 + 1 + {n — l)/2 — 6 = 1/2 — 6 for every (5 > 0. Iterated regularity 
with respect to this space then follows by interpolation. □ 

We now use this special case to prove Theorem 4.6. 

Proof of Theorem 4- 6- We begin by sharpening the regularity results of Theorem I.l 
to include uniformity and regularity in the location of the pole. 

Consider the Schwartz kernel of the fundamental solution, E{t,x,y;x',y'). For 
fixed {x',y') e (0,e) x dX for e > small enough and small t and x this is the 
distribution discussed in Theorem I.l. Moreover, by uniqueness of the solution to 
the Cauchy problem it depends continuously on {x',y'). It follows that the results 
described there hold uniformly in {x',y'). In particular near some fixed point t = 
X + x' with y, y' such that it is not on the direct front. 




5fn = Uit)iO,i5rn) 



(0, 



i6rn) e if-"/^+^-'(X°) ® H-''/^-%X°) 



{{t - x')dt + xd^fD'^E e iL, V k, a. 
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Here we have given up approximately half a derivative in the x,y variables, and 

continuity in x' , y' to settle for iterative regularity with respect to in all variables. 
By the symmetry and t-translation invariance of the problem E[t,x,y;x',y') = 
E{—t, x', y'; X, y) so it also follows that 

{{t - x)dt + x'd^.fD'^E e lL, V fc, a 

in the same region. All these vector fields commute, so by interpolation we deduce 
regularity with respect to all the vector fields simultaneously. Now, together with 
the wave operator itself, in both sets of variables, these symbols of these operators 
define the (two components of) the conormal bundle to {t = x + x'}. Thus we 
deduce that E is, away from the direct front and for small t, x, x' > 0, conormal 
with respect to this hyperplane. 

Although we have only shown iterative regularity in L^, iterative regularity in 

--e 

-^loc (R X X X) in this set follows by interpolation. Hence, a fortiori, E is 
conormal with respect to the diffracted front with iterate regularity in H^~'^(M. x 
X X X). 

Now, consider an admissible solution to the wave equation satisfying the hy- 
potheses of Theorem 4.6. Using Theorem 4.4 and a partition of unity, we may 
assume that the Cauchy data (u(0), D(m(0)) is identically zero in a neighborhood 
of (the projection to X of) all points geometrically related to p € R±. q{J.) and that 
at non geometrically related points it is supported in a microlocal neighborhood of 
R±,i(t). The nonfocusing condition then implies that there exists TV e N such that 
locally near x = t, 

u{t, X, y) 

I ^^^^^^^^^^uo{x',y')^dy' + I E{t,x,y,x\y')u,{x\y')^dy' 

we obtain the desired boundedness, by regarding E and dE/dt locally near the 
diffracted front as the kernels of Fourier integral operator of order e resp. 1 + e on 
M+, smoothly parametrized by t, with values in ^~°°{Y). □ 

The explicit construction of the fundamental solution in the product case by 
Cheeger- Taylor [2, 3] can be used to show that Theorem 4.6 cannot be strengthened 
by omitting the hypothesis of the nonfocusing condition. In particular, in Example 
4.1 of [2] the authors show that on X = R+ x the cone over the circle of 
circumference Air (with coordinate 9), the fundamental solution E{t, x, 0, x' , 9') has 
a jump discontinuity across the diffracted wavefront, where the value of the jump 
varies smoothly with 9,9', at least in \9 — 9'\ < tt. Now let (j){9) be a smooth 
function supported in \9\ < 5. The function 

/27r 
E{t,x,9,x',6')(l){9')d0' 
-27r 

is also an admissible solution to the wave equation, lying in C(M;r>i/2_e) for all 

e > 0. For t > x' , it has geometrically propagated singularities contained in [x = 
t - x',9 G {-S,d) ± tt). On the other hand, for all 9' G suppc/), E{t,x,9,x' ,9') 
has a jump discontinuity along x = t — x' , \9\ < it — 5, hence v also has a jump 
discontinuity this surface, i.e. this diffracted singularity is in 

for all e > 

0. Hence the diffracted singularity for v is no smoother than the geometrically 
propagated singularities. 



WAVE EQUATION ON CONIC MANIFOLDS 



53 



References 

[1] J. Briining and V.W. Guillcmin (Editors), Fourier integral operators, Springer- Verlag, Berlin, 

Heidelberg, New York, Tokyo, 1994. 
[2] Jeff Cheeger and Michael Taylor, On the diffraction of waves by conical singularities. I, 

Comm. Pure Appl. Math. 35 (1982), no. 3, 275-331, MR84h:35091a. 

[3] , On the diffraction of waves by conical singularities. II, Comm. Pure Appl. Math. 

35 (1982), no. 4, 487-529, MR84h:35091b. 
[4] J.J. Duistcrmaat and L. Hormandcr, Fourier integral operators, II, Acta Math. 128 (1972), 
183-269. 

[5] F. G. Fricdlandcr, Sound pulses, Cambridge University Press, New York, 1958. MR 20 #3703 
[6] Patrick Gerard and Gilles Lebeau, Diffusion d'une onde par un coin, J. Amer. Math. Soc. 6 
(1993), no. 2, 341-424. 

[7] M. W. Hirsch, C. C. Pugh, and M. Shub, Invariant manifolds. Springer- Verlag, Berlin, 1977, 

Lecture Notes in Mathematics, Vol. 583. MR 58 18595 
[8] L. Hormander, Fourier integral operators, I, Acta Math. 127 (1971), 79—183, See also [1]. 
[9] Lars Hormandcr, On the existence and the regularity of solutions of linear pseudo-differential 
equations, Enseigncmcnt Math. (2) 17 (1971), 99-163. MR 48 #9458 
[10] D. lagolnitzer. Appendix: Microlocal essential support of a distribution and decomposition 
theorems — an introduction, Hyperfunctions and theoretical physics (Rencontre, Nice, 1973; 
dedie a la memoire de A. Martineau), Springer, Berlin, 1975, pp. 121—132. Lecture Notes in 
Math., Vol. 449. MR 52 11583 
[11] M. Kalka and A. Menikoff, The wave equation on a cone. Comm. Partial Differential Equa- 
tions 7 (1982), no. 3, 223-278, MR83j:58110. 
[12] G. Lebeau, Propagation des ondes dans les varietes a coins, Seminaire sur les Equations aux 
Dcrivccs Partielles, 1995-1996, Ecole Polytech., Palaiseau, 1996, MR98m:58137, pp. Exp. No. 
XVI, 20. 

[13] Gillos Lebeau, Propagation des ondes dans les varietes d coins, Ann. Sci. Ecole Norm. Sup. 

(4) 30 (1997), no. 4, 429-497, MR98d:58183. 
[14] R. Mazzeo, Elliptic theory of differential edge operators I, Comm. in P.D.E. 16 (1991), 1615- 

1664. 

[15] R. B. Melrose and J. Wunsch, Singularities and the wave equation on conic spaces. National 
Research Symposium on Geometric Analysis and Applications (A. Isaev et al., ed.), vol. 39, 
Centre for Mathematics and its Applications, Australian National University, 2001. 

[16] R.B. Melrose, Spectral and scattering theory for the Laplacian on asymptotically Euclidian 
spaces, Spectral and scattering theory (Sanda, 1992) (M. Ikawa, cd.). Marcel Dckkcr, 1994, 
pp. 85-130. 

[17] Richard B. Melrose, The Atiy ah- Patodi- Singer index theorem, A K Peters Ltd., Wellesley, 

MA, 1993. MR 96g:58180 
[18] Detlef Miiller and Andreas Seeger, Inequalities for spherically symmetric solutions of the 

wave equation. Math. Z. 218 (1995), no. 3, 417-426. MR 96e:35093 
[19] Michel Rouleux, Diffraction analytique sur une variete a singularity conique, Comm. Partial 

Differential Equations 11 (1986), no. 9, 947-988. MR 88g:58182 
[20] B.-W. Schulzc, Boundary value problems and edge pseudo- differential operators, Microlocal 

analysis and spectral theory (Lucca, 1996), Kluwer Acad. Publ., Dordrecht, 1997, pp. 165- 

226. MR 98d:47112 

[21] Johannes Sjostrand, Singularites analytiques microlocales, Astcrisque, 95, Soc. Math. France, 

Paris, 1982, pp. 1-166. MR 84m:58151 
[22] A. Sommerfeld, Mathematische theorie der diffraktion. Math. Annalen 47 (1896), 317-374. 
[23] Michael E. Taylor, Pseudodifferential operators, Princeton University Press, Princeton, N.J., 

1981. MR 82i:35172 

[24] Jared Wunsch and Maciej Zworski, The FBI transform on compact C°° manifolds. Trans. 
Amer. Math. Soc. 353 (2001), no. 3, 1151-1167 (electronic). MR 1 804 416 

Department of Mathematics, Massachusetts Institute of Technology 
E-mail address: rbmamath . mlt . edu 

Department of Mathematics, SUNY at Stony Brook 
E-mail address: jwunschSmath. sunysb.edu 



